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EXECUTIVE  SUMMARY 


As  a result  of  industrial  development  in  the  Athabasca  oil 
sands  area  of  northeastern  Alberta,  concerns  have  been  raised 
regarding  potential  impacts  on  surface  water  quality  in  the  area.  It 
is  of  particular  concern  that  the  cumulative  effects  of  increasing 
development  of  the  oil  sands  may  create  hydrocarbon  and  metal  pollu- 
tion problems  in  the  Athabasca  River. 

With  the  aims  of  evaluating  baseline  water  quality  conditions 
and  identifying  areas  that  may  have  been  affected  by  existing  indus- 
trial developments,  water  quality  of  the  Athabasca  River  within  the 
oil  sands  area  has  been  monitored  since  1976.  In  addition,  several 
studies  of  benthic  invertebrate  communities  have  been  conducted  and 
these  provided  primarily  descriptive  information. 

This  report  presents  a detailed  statistical  analysis  of  water 
quality  and  benthic  invertebrate  data  from  previous  studies  on  the 
Athabasca  River.  Benthic  invertebrate  data  were  from  a study 
conducted  in  1981  on  the  Athabasca  River  between  Fort  McMurray  and 
the  Tar  River  confluence.  Water  quality  data  for  the  same  area  were 
obtained  from  the  NAQUADAT  water  quality  data  base  for  the  period 
1976  to  1983.  The  study  area  includes  a 75  km  section  of  the  Atha- 
basca River  extending  from  the  confluence  with  the  Horse  River 
upstream  of  Fort  McMurray,  downstream  to  the  confluence  with  the  Tar 
River,  approximately  40  km  downstream  from  the  Suncor  extraction  and 
upgrading  plant.  Six  water  quality  monitoring  stations  and  eight 
benthic  invertebrate  sampling  stations  were  located  within  the  study 
area. 

The  relationships  among  various  water  quality  parameters  were 
examined  using  the  methods  of  principal  component  analysis.  Princi- 
pal components  were  also  used  to  describe  associations  of  benthic 
invertebrate  taxa  and  for  transformation  of  abundance  data  prior  to 
making  statistical  comparisons  among  sampling  stations.  The 
relationships  between  water  quality  and  benthic  invertebrates  were 
examined  by  determining  the  correlations  of  benthic  invertebrate 
principal  components  with  water  quality  principal  components. 
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Strong  correlations  among  several  of  the  water  quality  para- 
meters were  apparent.  Most  of  the  major  ions  (calcium,  magnesium, 
sodium,  chloride,  sulphate,  and  bicarbonate)  were  strongly  correlated 
with  each  other  and  with  specific  conductance,  total  alkalinity,  and 
filterable  residue.  Potassium  concentrations  were  independent  of  the 
concentrations  of  other  major  ions.  There  were  also  components  of 
sodium  and  chloride,  attributable  to  the  clerwater  River  and  other 
east  bank  tributaries,  that  were  not  strongly  related  to  the  concen- 
trations of  other  ions.  All  of  the  metals  except  lead  and  mercury 
were  associated  with  non-f ilterable  residue  and  total  phosphate. 

With  respect  to  those  water  quality  parameters  included  in  the 
analyses,  there  was  no  evidence  that  effluent  from  the  Suncor  plant 
had  a large  or  consistent  effect  on  the  water  quality  of  the  Atha- 
basca River  downstream  from  the  development.  The  suggestion  is  that 
any  effects  on  water  quality  were  short-term  in  nature  and  did  not 
result  in  changes  that  persisted  for  long  periods  of  time. 

The  major  differences  in  water  quality  and  in  benthic  inverte- 
brate abundance  and  community  composition  within  the  study  area  were 
between  the  left  and  right  sides  of  the  Athabasca  River.  These 
differences  were  considered  to  be  due  primarily  to  the  influences  of 
the  Clearwater  River  and  other  east  bank  tributaries.  Some  observed 
differences  in  benthic  invertebrate  abundance  and  community  composi- 
tion within  the  study  area  may  be  related  to  nutrient  enrichment  from 
the  Fort  McMurray  sewage  effluent.  There  was  no  evidence  of  large 
differences  in  benthic  invertebrate  populations  between  stations 
immediately  upstream  and  downstream  of  the  Suncor  development  that 
could  be  attributed  to  the  Suncor  effluent. 

The  abundance  of  a variety  of  benthic  invertebrate  taxa 
appeared  to  be  correlated  with  several  water  quality  parameters.  In 
most  cases,  these  correlations  were  with  water  quality  parameters 
that  show  differences  between  the  left  and  right  sides  of  the  Atha- 
basca River.  The  observed  correlations  therefore  may  be  due  to  drif- 
ting of,  and  colonization  by,  invertebrates  from  the  Clearwater 


River,  and  possibly  other  right  bank  tributaries,  rather  than  to  any 
direct  influence  of  water  quality. 

The  results  of  statistical  analyses  of  the  benthic  inverte- 
brate data  presented  in  this  report  must  be  considered  tentative  or 
inconclusive.  A non-random  selection  procedure  was  used  during 
collection  of  the  benthic  invertebrate  samples  and  this  adversely 
affects  the  validity  of  statistical  analyses.  The  reliability  of  the 
results  is  therefore  questionable. 
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ABSTRACT 

A detailed  statistical  analysis  of  water  quality  and  benthic 
invertebrate  data  from  previous  studies  on  the  Athabasca  River  was 
undertaken.  Benthic  invertebrate  data  were  from  a study  conducted  in 
1981  on  the  Athabasca  River  between  Fort  McMurray  and  the  Tar  River 
confluence.  Water  quality  data  for  the  period  1976  to  1983  were 
obtained  from  the  NAQUADAT  water  quality  data  base. 

Principal  component  analysis  was  used  to  examine  the  relation- 
ships among  various  water  quality  parameters  and  was  also  used  for 
transformation  of  benthic  invertebrate  abundance  data  prior  to  making 
statistical  comparisons  among  sampling  stations.  The  relationships 
between  water  quality  and  benthic  invertebrates  were  examined  using 
correlation  analyses. 

Strong  correlations  among  several  of  the  water  quality  para- 
meters were  apparent.  Most  of  the  major  ions  were  correlated  with 
each  other  and  with  specific  conductance  and  filterable  residue. 
Potassium  concentrations  were  independent  of  the  concentrations  of 
other  major  ions.  All  of  the  metals  except  lead  and  mercury  were 
associated  with  non-f ilterable  residue  and  total  phosphate. 

The  major  differences  in  water  quality  and  in  benthic  inverte- 
brate abundance  and  community  composition  within  the  study  area  were 
between  the  left  and  right  sides  of  the  Athabasca  River.  These 
differences  were  considered  to  be  due  primarily  to  the  influences  of 
the  Clearwater  River  and  other  east  bank  tributaries.  Due  to  the  use 
of  a non-random  selection  procedure  during  collection  of  benthic 
invertebrate  samples,  the  results  of  statistical  analyses  based  on 
these  data  must  be  considered  inconclusive. 
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INTRODUCTION 


As  a result  of  industrial  development  in  the  Athabasca  oil 
sands  area  of  northeastern  Alberta,  concerns  have  been  raised 
regarding  potential  impacts  on  surface  water  quality  in  the  area.  It 
is  of  particular  concern  that  the  cumulative  effects  of  increasing 
development  of  the  oil  sands  may  create  hydrocarbon  and  metal 
pollution  problems  in  the  Athabasca  River. 

With  the  aims  of  evaluating  baseline  water  quality  condi- 
tions and  identifying  areas  that  may  have  been  affected  by  existing 
industrial  developments,  water  quality  of  the  Athabasca  River  within 
the  oil  sands  area  has  been  monitored  since  1976.  Water  quality  data 
are  maintained  on  the  NAQUADAT  storage  and  retrieval  system,  and  the 
results  of  the  monitoring  program  up  to  1981  have  been  reported  by 
Akena  (1982). 

In  addition,  studies  of  benthic  invertebrate  communities 
have  been  conducted  at  infrequent  intervals  (Flannagan  1976;  McCart 
et  al.  1977;  Barton  and  Lock  1979;  Barton  1980;  Boerger  1983).  These 
studies  provided  primarily  descriptive  information,  and  there  has 
been  no  attempt  to  examine  relationships  between  water  quality  and 
benthic  invertebrate  abundance  and  community  composition. 

The  present  study  is  concerned  with  a 75  km  section  of  the 
Athabasca  River  extending  from  the  confluence  with  the  Horse  River 
upstream  of  Fort  McMurray,  downstream  to  the  confluence  with  the  Tar 
River,  approximately  40  km  downstream  from  the  Suncor  extraction  and 
upgrading  plant.  Wastewater  effluents  entering  the  river  in  this 
area  include  effluent  from  the  Fort  McMurray  sewage  treatement  plant 
located  on  the  west  side  of  the  river,  the  Suncor  plant  discharges, 
which  also  enter  on  the  west  side  of  the  river  across  from  the  mouth 
of  the  Steepbank  River,  and  water  from  the  Syncrude  diversion,  which 
enters  the  Athabasca  River  through  Poplar  Creek,  approximately  9 km 
upstream  of  the  Suncor  effluent  outfall.  Water  from  the  Syncrude 
diversion  is  reported  to  be  a minor  influence,  but  may  contribute 
some  salt  loading.  Effluent  from  the  sewage  treatment  plant  contri- 
butes significant  nutrient  loadings  (orthophosphate  phosphorus  and 
ammonia)  as  well  as  a variety  of  organic  compounds  including  oil. 
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grease,  phenolic  material,  and  other  aliphatic  and  aromatic  hydro- 
carbons to  the  Athabasca  River  (Akena  1982).  The  Suncor  plant  dis- 
charges ammonia,  aliphatic  and  aromatic  hydrocarbons,  oils,  greases, 
phenol,  and  sulphide  into  the  river  (Akena  1982;  Boerger  1983).  In 
addition,  hydrocarbons  enter  the  river  from  natural  sources  as  a 
result  of  weathering  of  oil  sands  that  are  exposed  at  many  locations 
along  the  banks  of  the  Athabasca  River  and  the  banks  of  its  tribu- 
taries. Discharge  from  the  sewage  lagoons  at  Fort  McMurray  occurs 
every  month  throughout  the  year.  Release  from  sewage  lagoons  and 
Suncor  occurs  only  during  the  open  water  season,  and  the  releases  may 
be  interrupted. 

The  major  tributary  in  the  study  area  is  the  Clearwater  River, 
which  enters  the  Athabasca  River  at  Fort  McMurray.  The  Clearwater 
River  contributes  13  to  29%  of  the  annual  flow  of  the  Athabasca 
River,  and  no  other  tributary  in  the  study  area  accounts  for  more 
than  30%  of  the  mean  annual  flow  (Boerger  1983).  There  is  poor 
lateral  mixing  of  the  Athabasca  River  downstream  from  Fort  McMurray, 
and  water  from  the  Clearwater  River  remains  on  the  east  side  as  far 
downstream  as  the  Suncor  plant,  where  there  is  a sharp  bend  in  the 
river  (Boerger  1983).  Marked  differences  in  sodium  and  chloride  con- 
centrations between  the  east  and  west  sides  of  the  Athabasca  River 
have  been  observed  downstream  from  Fort  McMurray  and  were  attributed 
to  higher  sodium  chloride  levels  in  east  bank  tributary  streams  (par- 
ticularly the  Clearwater  River)  and  to  inflows  from  saline  reservoirs 
located  along  the  east  bank  of  the  Athabasca  River  (Akena  1982). 

This  report  presents  an  analysis  of  Athabasca  River  water 
quality  data  from  1976  to  1983  and  of  benthic  invertebrate  data 
collected  during  1981.  Correlations  among  selected  water  quality 
parameters  are  examined  and  seasonal  and  spatial  differences  in  water 
quality  are  described.  Benthic  invertebrate  abundance  and  distribu- 
tion data  (from  Boerger  1983)  are  analysed  in  detail  with  reference 
to  water  quality  and  to  sources  of  wastewater  effluents.  Correla- 
tions between  water  quality  and  benthic  invertebrate  abundance  and 
community  composition  are  examined. 
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2.  MATERIALS  AND  METHODS 

2.1  DATA  SOURCES 

Water  quality  data  for  seven  stations  on  the  Athabasca 
River  were  obtained  from  the  NAQUADAT  storage  and  retrieval  system 
for  the  period  1976  to  1983.  These  data  were  provided  on  magnetic 
tape  by  the  Water  Quality  Control  Branch,  Alberta  Environment. 

Benthic  invertebrate  sampling  was  conducted  at  eight 
stations  in  the  study  area  at  two-week  intervals  during  the  spring 
and  summer  of  1981.  The  results  of  that  study  have  been  reported  by 
Boerger  (1983).  Benthic  invertebrate  data  used  for  the  present  study 
were  taken  from  an  appendix  to  that  report. 

2.2  SAMPLING  LOCATIONS  AND  SCHEDULES 
2.2.1  Water  Quality 

Six  stations  on  the  Athabasca  River  within  the  study  area 
have  been  sampled  as  part  of  the  regular  water  quality  monitoring 
program  conducted  by  the  Pollution  Control  Division,  Alberta  Environ- 
ment. The  station  identification  numbers  and  their  locations  as 
given  by  Akena  (1982)  are  as  follows:  Station  CC0012,  100  m upstream 
from  the  confluence  with  the  Horse  River;  Station  DA0203,  at  Mile  19 
(just  above  the  confluence  with  Poplar  Creek);  Station  DA0205,  at 
Mile  29.8  (approximately  7 km  downstream  from  the  Suncor  effluent 
outfall);  Station  DA0206,  at  Mile  34.5  (just  upstream  from  the 
confluence  with  the  Muskeg  River);  Station  DA0207,  at  Fort  MacKay 
(below  the  confluence  with  the  MacKay  River);  Station  DA0208,  at  Mile 
52.4  (approximately  4 km  downstream  from  the  confluence  with  the  Tar 
River).  One  additional  station  (DD0010,  at  the  WSC  gauge  at  Embarras 
Airport),  which  is  located  a considerable  distance  downstream  from 
the  study  area,  was  also  included  in  analyses  examining  the  correla- 
tions among  water  quality  parameters. 

The  locations  of  three  of  the  stations  (DA0205,  DA0207,  and 
DA0208)  apparently  have  been  changed  from  those  given  above.  The 
following  information  on  the  locations  at  which  these  stations  are 
currently  sampled  was  provided  by  Lynda  Corkum  (Pollution  Control 
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Division,  Water  Quality  Control  Branch,  Alberta  Environment). 
Station  DA0205  is  now  located  immediately  below  the  Suncor  effluent 
outfall,  approximately  7 km  upstream  from  its  former  location. 
Station  DA0207  is  currently  sampled  just  above  the  confluence  with 
the  MacKay  River  rather  than  below  it.  The  present  location  of 
Station  DA0208  is  upstream  from  the  confluence  with  the  Tar  River,  4 
to  5 km  upstream  from  the  earlier  location.  We  were  not  able  to 
ascertain  when  these  changes  were  made.  However,  statements  made  by 
Boerger  (1983)  regarding  the  correspondence  between  his  benthic 
invertebrate  sampling  stations  and  water  quality  stations  suggest 
that  the  locations  of  stations  DA0207  and  DA0208  in  1981  were  the 
same  as  the  present  locations.  The  locations  of  the  water  quality 
monitoring  sites  are  shown  on  Figure  1. 

All  of  the  water  quality  stations  except  DD0010  were 
sampled  regularly  during  1976  and  into  early  1977.  Stations  CC0012 
and  DA0207  were  also  monitored  from  mid-1977  through  1983.  Station 
DD0010  was  monitored  regularly  from  late  1978  through  1983.  Station 
DA0206  was  not  sampled  from  1977  through  1979,  but  has  been  sampled 
regularly  since  1980.  Stations  DA0203  and  DA0208  were  not  monitored 
from  1977  through  1980.  Station  DA0203  was  sampled  regularly  from 
1981  through  1983  and  Station  DA0208  was  sampled  regularly  in  1981 
and  1983,  but  not  for  most  of  1982.  Station  DA0205  was  not  monitored 
from  1977  through  1981,  but  was  monitored  regularly  during  1982  and 
1983.  During  the  periods  in  which  a station  was  monitored,  sampling 
was  generally  on  a monthly  basis,  but  there  are  often  gaps  of  several 
months  (usually  during  winter)  in  the  records  for  stations  CC0012  and 
DD0010.  In  the  spring  and  summer  of  1981,  stations  CC0012,  DA0203, 
DA0206,  DA0207,  and  DA0208  were  sampled  at  two-week  intervals  on 
dates  coinciding  with  the  benthic  invertebrate  sampling  program. 

During  the  period  1976  to  1983,  538  water  samples  were 
taken  from  the  stations  and  68  water  quality  parameters  were 
measured.  Not  all  parameters  were  determined  for  every  sample, 
however,  and  for  any  one  sample,  the  number  of  parameters  measured 
was  usually  much  less. 
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Figure  1.  Map  of  the  study  area  showing  locations  of  water  quality 
monitoring  sites  and  benthic  invertebrate  sampling 
stations . 
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After  1983,  the  monitoring  program  at  most  of  the  stations 
generally  included  samples  taken  on  both  the  left  and  right  sides  of 
the  river  (looking  downstream).  Stations  CC0012  and  DD0010  were 
usually  sampled  only  at  mid-channel.  The  location  in  the  river 
channel  was  indicated  by  distances  to  the  left  and  right  banks,  and 
the  distance  to  the  closest  bank  for  samples  apparently  taken  as  left 
or  right  bank  samples  ranged  from  1 m to  approximately  one-third  of 
the  river  width.  In  the  present  study,  samples  located  anywhere 
within  the  left  34%  of  the  river  width  were  considered  left  bank 
samples,  those  within  the  right  34%  were  considered  right  bank 
samples,  and  those  within  the  central  32%  were  considered  mid-channel 
samples.  For  105  of  the  538  water  samples,  there  was  insufficient 
information  in  the  NAQUADAT  files  to  determine  sampling  location. 
These  samples  were  assumed  to  be  mid-channel  samples. 

2.2.2  Benthic  Invertebrates 

The  eight  stations  on  the  Athabasca  River  sampled  for 
benthic  invertebrates  by  Boerger  (1983)  during  1981  were  as  follows: 
Station  1,  between  the  Tar  River  and  Ells  River  confluences;  Station 
2,  above  the  confluence  with  the  MacKay  River;  Station  3,  approxi- 
mately 2 km  upstream  from  the  confluence  with  the  Muskeg  River; 
Station  4,  approximately  4 km  downstream  from  the  Suncor  effluent 
outfall;  Station  5,  upstream  of  the  Suncor  development  and  below  the 
confluence  with  Poplar  Creek;  Station  6,  approximately  5 km  upstream 
from  the  confluence  with  Poplar  Creek;  Station  7,  approximately  6 km 
downstream  from  Fort  McMurray;  Station  8,  near  the  confluence  with 
the  Horse  River,  upstream  of  Fort  McMurray.  The  locations  of  these 
stations  are  indicated  on  Figure  1.  At  each  station,  samples  were 
taken  from  gravel  bars  on  both  the  east  (right)  and  west  (left)  sides 
of  the  river,  giving  a total  of  16  sites  (designated  IE  and  1W 
through  8E  and  8W) . 

A cylinder  sampler  was  used,  and  three  samples  were  taken 
from  each  of  the  16  sites  at  each  sampling  interval.  Sampling  was 
conducted  at  two-week  intervals  on  the  following  dates  in  1981:  May 
13  and  14,  May  28  and  29,  June  9 and  10,  June  23  and  24,  Duly  7 and 
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8,  July  21  and  22,  August  and  August  18  and  19.  The  August 
sampling  was  incomplete,  and  the  data  for  that  interval  were  not  used 
in  the  present  study. 

Five  of  the  stations  were  located  at  or  near  water  quality 
monitoring  stations.  Benthic  invertebrate  stations  1,  2,  3,  6,  and  8 
correspond  to  water  quality  stations  DA0208,  DA0207,  DA0206,  DA0203, 
and  CC0012,  respectively. 

2.3  STATISTICAL  METHODS 

2.3.1  Principal  Component  Analysis 

Principal  component  analysis  (PCA)  was  used  as  an  initial 
step  in  the  analysis  of  both  water  quality  and  benthic  invertebrate 
data.  Application  of  PCA  to  water  quality  data  was  used  to  examine 
the  relationships  among  various  water  quality  parameters  and  to 
reduce  the  number  of  variables  for  subsequent  analyses.  With  benthic 
invertebrate  data,  PCA  was  used  primarily  as  a data  reduction  trans- 
formation to  allow  use  of  fewer  variables  in  subsequent  statistical 
analyses.  All  principal  component  analyses  were  performed  using  the 
BMDP  Statistical  Software  package  of  computer  programs  (Dixon  et  al. 
1981). 

The  method  of  principal  component  analysis  is  described  in 
detail  by  Morrison  (1976)  and  Gohnson  and  Wichern  (1982).  Green 
(1979)  discusses  applications  of  PCA  in  environmental  studies.  The 
following  introductory  description  is  adapted  from  Marriott  (197A-) 
and  Morrison  (1976). 

Principal  component  analysis  transforms  a set  of  variables 
X^...Xp  to  a new  set  of  variables  Y'|...Yp  having  the  following 
properties: 

1.  Each  new  variable  (principal  component)  is  a linear 
combination  of  the  original  variables;  for  example 

Yi  = aiiXi  + a^2^2  + •••  + ^ip^p 

2.  The  sum  of  squares  of  the  component  coefficients 
j=1  ...  p,  is  unity . 
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3.  Of  all  possible  combinations  of  this  type,  (the 
first  principal  component)  has  the  greatest  variance 
(i.e.  accounts  for  the  greatest  proportion  of  the  total 
variance) . 

4.  Of  all  possible  combinations  of  this  type  not  correla- 
ted with  Yi,  Y2  has  the  greatest  variance;  Y3  has  the 
greatest  variance  of  combinations  not  correlated  with 
Y-j  and  Y£,  and  so  on  until  the  complete  set  of  princi- 
pal components,  from  Y^  to  Yp,  has  been  defined. 

The  principal  components  so  defined  are  not  correlated  with 
each  other  and  are  arranged  in  order  of  decreasing  variance 
explained.  Frequently,  the  first  few  principal  components  will 
account  for  a large  proportion  of  the  variability  in  the  original 
data,  and  these  few  can  then  be  used  in  further  analyses,  in  place  of 
the  original  variables,  with  minimal  loss  of  information. 

Interpretation  of  principal  components  is  based  on  the 
loadings  of  each  of  the  original  variables  on  each  of  the  principal 
components.  These  loadings  are  the  correlations  of  the  original 
variables  with  the  principal  components,  which  are  generally  more 
useful  than  the  component  coefficients  for  interpretive  purposes 
(Johnson  and  Wichern  1982).  Variables  with  high  loadings  on  the  same 
principal  component  tend  to  be  highly  correlated  with  each  other,  and 
each  principal  component  is  interpreted  according  to  the  magnitudes 
of  the  loadings  associated  with  it.  In  some  cases,  the  principal 

components  clearly  may  represent  some  identifiable  factors  such  as 
various  environmental  influences.  Although  this  kind  of 

interpretation  usually  is  not  possible  for  benthic  invertebrates, 
each  principal  component  of  benthic  invertebrate  abundance  data  still 
may  still  be  interpreted  in  terms  of  the  abundances  of  the 
constitutent  taxa. 

For  the  purpose  of  interpreting  principal  components  when 
the  number  of  original  variables  is  large,  it  is  useful  to  consider 
only  those  variables  with  relatively  large  (absolute  value)  loadings 
on  each  principal  component.  Determination  of  what  represents  a 
large  loading  is  arbitrary  and  will  depend  to  some  extent  on  the 
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number  of  original  variables  and  the  sample  size.  When  interpreting 
principal  components  in  the  present  study,  we  have  considered  only 
loadings  greater  than  0.25  absolute  value.  This  approach  is  used 
only  for  convenience  of  interpretation.  Since  principal  component 
scores  are  calculated  using  all  of  the  original  variables,  subsequent 
analyses  based  on  these  scores  are  not  affected  by  the  criteria  used 
to  select  large  loadings. 

2. 3. 1.1  Water  quality  PCA.  The  determination  of  which  water 
quality  variables  to  include  in  a principal  component  analysis  was 
based  on  several  considerations.  Since  variables  that  are  calculated 
from  other  variables  (e.g.,  total  hardness)  are  redundant,  these  were 
not  included.  In  addition,  variables  that  were  almost  always  at  or 
below  the  detection  limit  (e.g.,  carbonate)  provided  little 
information  and  these  were  also  omitted  from  the  analysis.  PCA 
requires  measurements  of  all  variables  for  all  samples.  Not  all 
variables  were  determined  for  every  water  sample,  however,  and  some 
variables  were  measured  infrequently.  The  selection  of  variables  was 
necessarily,  therefore,  a compromise  between  the  desire  to  include  as 
many  variables  as  possible  and  the  desire  to  include  as  many  water 
samples  as  possible. 

After  examining  the  number  of  samples  included  when  several 
different  subsets  of  water  quality  variables  were  selected,  a set  of 
29  variables  was  chosen  as  representing  a reasonable  compromise  in 
terms  of  maximizing  both  the  number  of  variables  and  the  number  of 
samples.  These  29  variables  were  subjected  to  a principal  component 
analysis  and  3A-3  water  samples  were  included.  This  selection  proced- 
ure eliminated  195  samples,  including  all  of  those  taken  prior  to 
1978. 

For  the  purpose  of  examining  the  relationships  between 
water  quality  and  benthic  invertebrates,  it  was  desirable  to  include 
all  water  samples  taken  during  1981  at  locations  where  benthic  inver- 
tebrates were  also  sampled.  This  required  further  reduction  of  the 
number  of  water  quality  variables  to  a set  of  20.  A second  PCA 
therefore  was  performed  using  these  20  variables.  This  analysis 
included  435  water  samples. 
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Preliminary  examination  of  the  distributions  of  values  for 
each  of  the  water  quality  variables  suggested  that  a logarithmic 
transformation  would  be  appropriate  (to  more  closely  approximate 
normality)  for  all  variables  except  pH.  A log-jQ  transformation  was 
therefore  applied  to  all  variables  except  pH  prior  to  performing  a 
PCA.  Values  that  were  recorded  as  being  less  than  a detection  limit 
were  changed  to  a value  slightly  less  than  the  detection  limit  before 
transforming  to  logarithms.  For  example,  if  the  detection  limit  was 
0.02,  a value  of  0.019  was  assigned  to  measurements  reported  as  less 
than  that  limit.  This  assignment  of  values  was  done  by  the  Water 
Quality  Control  Branch  before  the  data  files  were  made  available. 

Principal  components  were  derived  from  the  correlation 
matrix  of  water  quality  variables  and  an  orthogonal  rotation  of  the 
loadings  matrix  was  applied  using  the  Varimax  criterion  (Harman 
1976).  The  purpose  of  the  rotation  was  to  allow  a simpler  interpre- 
tation of  the  principal  components.  Components  that  accounted  for 
more  than  3%  of  the  total  variance  were  retained  for  subsequent 
analysis. 


2. 3. 1.2  Benthic  invertebrates  PCA.  Because  the  composition  of  ben- 
thic invertebrate  communities  changes  substantially  on  a seasonal 
basis,  it  is  preferable  to  apply  the  method  of  PCA  to  data  collected 
over  a relatively  short  period  of  time.  Accordingly,  four  separate 
analyses  were  performed:  one  for  each  of  the  months  of  May,  3une, 
July,  and  August.  Data  from  the  two  sampling  periods  in  each  of  May, 
June,  and  July  were  pooled  for  these  analyses.  In  August,  data  was 
available  for  only  one  sampling  period. 

Rare  taxa  usually  provide  little  information  and  can  some- 
times create  computational  difficulties  in  PCA.  Therefore,  for  each 
analysis,  only  those  taxa  comprising  at  least  10%  of  any  one  sample, 
or  found  in  at  least  10%  of  the-  samples,  were  included.  Analyses 
were  performed  using  taxon  abundances  (number/m^)  transformed  to 
log^iQ  (abundance  +1).  Principal  components  were  derived  from  the 
covariance  matrix  of  transformed  abundances.  This  is  the  preferred 
method  of  PCA  when  all  variables  are  measured  in  the  same  units 
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(Morrison  1976;  Green  1979).  Because  orthogonal  rotation  of  the 
loadings  matrix  did  not  yield  a simpler  interpretation  of  the  princi- 
pal components,  the  unrotated  components  were  used.  Principal 
components  that  accounted  for  more  than  A-%  of  the  total  variance  were 
retained  for  subsequent  analysis. 

2.3.2  Other  Statistical  Procedures 

Statistical  comparisons  of  sampling  stations,  with  respect 
to  benthic  invertebrates,  were  made  using  analysis  of  variance  and 
Student-Newman-Keuls  (S-N-K)  multiple  comparisons  procedures  (Sokal 
and  Rohlf  1969;  Winer  1971).  These  comparisons  were  based  on  the 
principal  components  of  taxon  abundances.  A one-way  analysis  of 
variance  was  first  used  to  test  each  principal  component  for  any 
differences  among  stations.  For  any  principal  component  that  showed 
a significant  difference  among  stations,  all  pairwise  comparisons  of 
stations  were  then  made  using  the  S-N-K  multiple  comparisons  method. 

The  relationships  between  water  quality  and  benthic  inver- 
tebrates were  examined  using  simple  correlation  analysis,  based  on 
correlations  among  principal  components.  The  Pearson  product-moment 
correlation  coefficient  (Sokal  and  Rohlf  1969)  was  computed  for  each 
benthic  invertebrate  principal  component  paired  with  each  water 
quality  principal  component,  and  tested  for  significance.  Examina- 
tion of  scattergrams  gave  no  indication  that  any  relationship  other 
than  linear  would  be  appropriate  for  any  pair  of  components. 

Correlation  analyses,  analysis  of  variance,  and  S-N-K 
multiple  comparisons  tests  all  were  performed  using  computer  programs 
of  the  Statistical  Package  for  the  Social  Sciences  (Nie  et  al.  1975). 


12 


3.  RESULTS  AND  DISCUSSION 

3.1  WATER  QUALITY 

3.1.1  Water  Quality,  1976  to  1983 

For  reasons  outlined  earlier  (Section  2. 3. 1.1),  it  was  not 
possible  to  include  all  68  water  quality  variables  in  a principal 
component  analysis.  The  data  for  all  variables  are,  however, 
presented  for  reference  purposes  in  Figures  16  to  83  (Appendix  7.1). 
There  is  one  figure  for  each  variable,  and  each  figure  includes  data 
for  the  seven  sampling  stations,  covering  the  period  1976  to  1983. 

The  results  of  principal  component  analysis  of  29  selected 
water  quality  variables  are  presented  in  Table  1,  which  contains  the 
principal  component  loadings  for  each  water  quality  variable.  Eight 
principal  components  (PCs)  accounted  for  78%  of  the  total  variance, 
PCI  accounted  for  36.^%  and  PC2  accounted  for  17.0%  of  the  total 
variance  while  the  amount  of  variance  explained  by  each  of  PCs 
3,  4,  5,  6,  7,  and  8 was  5.8%,  5.1%,  4.0%,  3.7%,  3.3%,  and  3.2%, 
respectively.  The  scores  for  a PC,  which  are  computed  for  each  water 
sample,  are  determined  primarily  by  the  variables  having  the  largest 
loadings  (absolute  value)  on  that  PC. 

The  first  principal  component  (PCI)  has  high  loadings  for 
specific  conductance,  total  alkalinity,  filterable  residue,  and  all 
the  major  ions  (although  the  loading  for  potassium  is  relatively 
low).  Nitrate  and  nitrite  nitrogen  and  ammonia  also  are  associated 
with  this  PC,  as  is  iron.  Non-f ilterable  residue  and  total  organic 
carbon  are  inversely  related  to  this  component  (indicated  by  the 
negative  signs  of  the  loadings).  PCI  clearly  represents  primarily 
dissolved  solids,  in  a general  way. 

PC2  represents  non-f ilterable  residue,  total  organic 
carbon,  total  phosphate,  and  all  the  metals  except  mercury.  The 
loading  for  lead  is  relatively  low,  however.  Ammonia  is  also  asso- 
ciated with  this  component  to  some  degree,  as  are  sodium  and 
chloride,  the  latter  two  in  a negative  manner.  A strong  association 
between  non-f ilterable  residue,  total  phosphate,  and  most  of  the 
metals  is  indicated. 


Table  1.  Principal  component  loadings  for  the  first  eight  principal  components  of  29  water  quality 
variables.  Loadings  greater  than  0.25  absolute  value  are  printed  in  boldface. 
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The  third  PC  is  primarily  a measure  of  sodium  and  chloride 
that  is  independent  of  PC1  (and  therefore  of  the  other  major  ions). 
Ammonia  and  nitrate  plus  nitrite  nitrogen  are  also  associated  with 
PC3  to  some  degree,  and  pH  is  inversely  related  to  it.  It  will 
become  clear  later  that  PC3  represents  primarily  the  higher  sodium 
and  chloride  levels  of  the  Clearwater  River  and,  possibly,  other 
right  bank  tributary  streams. 

Potassium  is  the  only  variable  with  a large  positive  load- 
ing on  PC4.  Vanadium  and  pH  are  inversely  related  to  this  PC.  The 
high  loading  for  potassium  on  PC4  indicates  that  the  concentration  of 
potassium  is  largely  independent  of  the  concentrations  of  the  other 
major  ions. 

PC5  has  high  loadings  for  oil  and  grease,  lead,  and  zinc, 
indicating  an  association  between  oil  and  grease  and  at  least  some 
forms  of  lead  and,  to  a lesser  extent,  zinc.  Cyanide  has  a strong 
negative  correlation  with  PC5,  and  this  is  the  only  PC  with  which 
cyanide  is  strongly  associated. 

Dissolved  orthophosphate  loads  heavily  on  PC6  and  total 
organic  carbon  and  lead  also  are  associated  with  this  component,  but 
to  a lower  degree.  PC6  therefore  represents  primarily  orthophosphate 
and  is  the  only  PC  with  which  orthophosphate  is  strongly  associated. 

PC7  represents  primarily  mercury,  and  the  loading  for  oil 
and  grease  on  this  PC  indicates  a negative  correlation  between 
mercury  and  oil  and  grease. 

Phenolic  material  and  oil  and  grease  are  the  only  variables 
associated  strongly  with  PCS.  This  component  is  therefore  primarily 
a measure  of  these  two  variables. 

Plots  of  the  eight  principal  components  are  presented  in 
Figures  2 to  9 for  all  seven  water  quality  monitoring  stations  during 
the  period  1978  to  1983.  All  samples  taken  prior  to  1978  were 
omitted  from  this  analysis  because  they  lacked  measurements  for  one 
or  more  of  the  water  quality  variables. 

As  would  be  expected,  PCI  (representing  dissolved  solids) 
shows  marked  seasonal  fluctuations  with  the  highest  values  occurring 
during  the  winter  months.  Differences  between  left  and  right  bank 
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samples  were  not  large  in  most  instances.  Values  of  PCI  were 
consistently  higher,  however,  for  left  bank  samples  at  Station  DA0205 
(below  the  Suncor  effluent  outfall)  in  1982  and  1983,  particularly 
during  winter.  Station  DA0203  also  had  higher  values  of  PCI  in  left 
bank  samples,  at  least  during  winter,  and  the  difference  may  there- 
fore reflect  only  differences  between  the  Athabasca  and  Clearwater 
rivers  rather  than  an  effect  of  the  Suncor  effluent.  In  general, 
there  were  no  substantial  differences  between  any  of  the  stations 
with  respect  to  PCI . 

The  highest  values  for  PC2,  which  represents  non-f ilterable 
residue,  total  organic  carbon,  total  phosphate,  and  several  metals, 
occurred  during  the  spring  and  summer  months  and  were  presumably 
related  to  high  discharges.  There  was  no  evidence  of  consistent 
differences  between  left  and  right  bank  samples  at  any  station,  and 
the  values  at  all  stations  were  usually  similar. 

PC3,  primarily  a measure  of  sodium  and  chloride,  showed 
much  higher  values  in  right  bank  samples  than  in  left  bank  samples  at 
both  stations  DA0203  and  DA0205,  and  clearly  reflects  the  influence 
of  the  Clearwater  River,  which  has  higher  sodium  and  chloride  concen- 
trations than  the  Athabasca  River  (Akena  1982).  At  stations  DA0206 
and  DA0207,  right  bank  samples  frequently  had  higher  PC3  values  than 
left  bank  samples,  but  the  differences  were  not  as  great  or  as 
consistent  as  at  stations  DA0203  and  DA0205.  The  differences 
observed  at  stations  DA0206  and  DA0207  may  be  the  result  of 
influences  of  right  bank  tributary  streams  other  than  the  Clearwater 
River,  or  of  saline  groundwater  sources.  Mixing  of  the  Athabasca 
River  is  reported  to  occur  at  a sharp  bend  in  the  river  of  the  Suncor 
dyke  (Boerger  1983).  The  large  differences  observed  at  Station 
DA0205  suggest,  however,  that  mixing  may  not  be  complete  at  this 
point,  which  is  a short  distance  downstream  from  the  dyke.  The 
highest  values  for  PC3  most  often  occurred  in  the  first  months  of  the 
year,  but  seasonal  fluctuations  were  not  as  regular  as  was  the  case 
with  PC1 . 

PC^,  which  is  related  primarily  to  potassium  (directly)  and 
vanadium  (inversely),  did  not  show  large  or  consistent  differences 
between  left  and  right  bank  samples  at  any  station,  and  values  did 


16 


1976  1977  1978  1979  1980  1981  1982  1983 

LEGEND 

Left  bank 

Water  Quality  right  bank 

Principal  Component  1 ml3”channel 


Figure  2.  Principal  component  scores  for  water  quality  PC1  at  seven 
stations  on  the  Athabasca  River. 


continued  . 
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Figure  3.  Principal*' component  scores  for  water  quality  PC2  at  seven 
stations  on  the  Athabasca  River. 


continued  . 
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Figure  3.  Concluded. 
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Figure  4.  Principal  component  scores  for  water  quality  PC3  at  seven 
stations  on  the  Athabasca  River. 
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Figure  Concluded. 
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Figure  5.  Principal  component  scores  for  water  quality  PC4  at  seven 
stations  on  the  Athabasca  River. 


continued  . 
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Figure  6.  Principal  component  scores  for  water  quality  PC5  at  seven 
stations  on  the  Athabasca  River. 
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Figure  7.  Principal  component  scores  for  water  quality  PC6  at  seven 
stations  on  the  Athabasca  River. 
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Figure  8.  Principal  component  scores  for  water  quality  PC7  at  seven 
stations  on  the  Athabasca  River. 
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Figure  9.  Principal  component  scores  for  water  quality  PCS  at  seven 
stations  on  the  Athabasca  River. 


continued  . 
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Figure  9.  Concluded. 
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not  appear  to  differ  markedly  between  stations.  Although  changes  in 
PC4  during  the  year  were  evident,  these  changes  were  not  consistently 
related  to  season  from  year  to  year. 

Values  for  PC5  (representing  lead,  oil  and  grease,  and 
zinc  on  the  positive  side,  and  cyanide  on  the  negative  side)  were 
generally  similar  at  all  stations,  and  there  was  no  clear  indication 
of  any  consistent  differences  between  left  and  right  bank  samples. 
Fluctuations  of  PC5  did  not  appear  to  be  seasonal  in  nature. 

Although  the  highest  values  for  PC6,  which  represents  prim- 
arily dissolved  orthophosphate  and,  to  a lesser  extent,  lead  and 
total  organic  carbon,  tended  to  occur  in  early  spring,  the  pattern  of 
seasonal  changes  was  not  very  regular.  There  was  no  evidence  of 
differences  between  stations,  and  PC6  was  usually  similar  for  left 
and  right  bank  samples  at  all  stations. 

PC7,  which  represents  primarily  mercury,  did  not  show 
regular  seasonal  changes  and,  except  for  occasional  peaks,  values 
were  similar  at  all  stations  and  in  left  and  right  bank  samples.  Two 
of  the  largest  values  for  PC7  occurred  at  Station  DA0203;  one  in  a 
right  bank  sample  in  mid-1982,  and  one  in  a left  bank  sample  in  May 
of  1983.  A similar  peak  was  observed  in  a left  bank  sample  at 
Station  DA0208  in  May  1983. 

Fluctuations  of  PC8  (representing  primarily  phenolic 

material  and  oil  and  grease)  did  not  show  a con^sistent  seasonal  pat- 

■ 

tern.  Values  for  left  and  right  bank  samples  were  usually  similar, 
and  no  consistent  differences  between  stations  were  apparent. 

3.1.2  Water  Quality,  1981 

In  order  to  assess  relationships  between  water  quality  and 
benthic  invertebrates,  it  was  desirable  to  include  in  analyses  all 
water  samples  taken  in  1981  at  locations  and  times  corresponding  to 
the  benthic  invertebrate  sampling  program.  For  this  purpose,  it  was 
necessary  to  restrict  the  water  quality  variables  considered  to  a 
group  of  20.  A second  principal  component  analysis  therefore  was 
performed  using  these  20  variables. 
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Results  of  a principal  component  analysis  of  the  selected 
variables  are  presented  in  Table  2,  which  contains  the  principal  com- 
ponent loadings  for  each  water  quality  variable.  The  PC  scores 
for  each  sample  are  determined  primarily  by  the  variables  having  the 
largest  (i.e.,  at  least  0.25  absolute  value)  loadings  on  each  PC. 
Six  principal  components  accounted  for  77%  of  the  total  variance  and, 
for  the  most  part,  results  were  consistent  with  the  results  of  the 
principal  component  analysis  of  29  water  quality  variables  presented 
in  Section  3.1.1.  PCI  accounted  for  39.5%  and  PC2  accounted  for 
16.6%  of  the  total  variance.  The  percentage  of  the  variance 
explained  by  each  of  PC3,  PC^,  PC5,  and  PC6  was  6.3%,  5.7%,  4.6%,  and 
4.3%,  respectively. 

The  first  principal  component  (PCI ) has  high  loadings  for 
specific  conductance,  total  alkalinity,  filterable  residue,  calcium, 
magnesium,  sodium,  and  chloride.  Nitrate  plus  nitrite  nitrogen  is 
also  associated  with  this  PC,  and  non-f ilterable  residue  and  total 
organic  carbon  are  negatively  correlated  with  it.  In  general,  PCI 
can  be  considered  to  represent  primarily  dissolved  solids. 

PC2  represents  non-f ilterable  residue,  total  phosphate, 
total  organic  carbon,  and  several  metals  (manganese,  copper,  zinc, 
and  vanadium) . Sodium  and  chloride  have  an  inverse  relationship  with 
this  PC.  The  loadings  for  PC2  indicate  a strong  relationship  between 
non-f ilterable  residue,  total  phosphate,  and  most  of  the  metals. 

The  third  PC  represents  primarily  potassium,  copper,  zinc, 
and  lead.  The  large  loading  for  potassium  on  PC3  and  the  relatively 
low  loading  for  potassium  on  PCI  indicate  that  the  concentration  of 
potassium  is  largely  independent  of  the  concentrations  of  the  other 
major  ions. 

PC4  has  high  loadings  for  sodium,  chloride,  ammonia,  and  to 
a lesser  extent,  nitrate*  plus  nitrite  nitrogen.  This  component  is 
primarily  a measure  of  sodium  and  chloride  (that  is  independent  of 
PCI)  and  ammonia.  As  was  the  case  with  the  third  principal  component 
in  the  analysis  of  Section  3.1.1,  PC4  reflects  the  influence  of  the 
Clearwater  River  and,  possibly,  other  right  bank  tributary  streams. 


Table  2.  Principal  component  loadings  for  the  first  six  principal  components  of  20  water  quality 
variables.  Loadings  greater  than  0.25  absolute  value  are  indicated  in  boldface. 
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The  variable  loading  most  heavily  on  PC5  is  dissolved 
orthophosphate,  and  this  is  the  only  PC  with  which  orthophosphate  is 
strongly  associated.  Total  organic  carbon  and  potassium  also  are 
associated  with  PC5  to  some  degree,  but  their  loadings  are  much 
smaller  than  the  loading  for  orthophosphate.  PC5  may  therefore  be 
considered  to  be  primarily  a measure  of  dissolved  orthophosphate. 

Phenolic  material  has  a very  high  loading  on  PC6,  and  the 
loadings  for  all  other  variables  on  this  component  are  very  low.  PC6 
therefore  essentially  represents  only  phenolic  material. 

Plots  of  the  six  principal  components  for  water  samples 
taken  during  1981  are  presented  in  Figures  10  to  15.  The  PC  scores 
presented  are  for  water  quality  monitoring  stations  CC0012,  DA0203, 
DA0206,  DA0207,  and  DA0208,  which  correspond  to  benthic  invertebrate 
sampling  stations  8,  6,  3,  2,  and  1,  respectively. 

Values  of  PCI,  which  represents  primarily  dissolved  solids, 
declined  from  highs  in  February  and  March  to  lows  during  the  spring 
and  summer  months,  and  then  began  increasing  from  September  on.  In 
general,  PCI  did  not  differ  markedly  between  stations,  but  some  dif- 
ferences between  left  and  right  bank  samples  were  apparent.  At 
Station  DA0203,  PCI  scores  were  higher  for  left  bank  samples  from  May 
until  late  August,  with  the  exception  of  one  sample  taken  in  late 
3une.  Similar  differences  were  also  evident  at  stations  DA0206  and 
DA0207  until  late  August,  again  with  the  exception  of  one  sampling 
date  in  Dune.  These  differences  probably  reflect  the  influence  of 
the  Clearwater  River  and  other  right  bank  tributaries.  They  do  not 
appear  to  be  related  to  sources  of  effluents. 

The  highest  values  for  PC2  (representing  primarily  non-fil- 
terable  residue,  total  phosphate,  and  several  metals)  occurred  during 
the  spring  and  summer  months  and  were  undoubtedly  associated  with 
periods  of  high  discharge.  In  general,  differences  between  stations 
were  not  large  or  consistent.  There  were  some  differences  between 
left  and  right  bank  samples,  however.  Except  on  one  date  in  3une, 
PC2  scores  were  consistently  higher  for  left  bank  samples  at  Station 
DA0203.  At  Station  DA0206,  PC2  was  higher  for  left  bank  samples  on 
all  dates  except  during  the  period  from  late  August  to  early 
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Figure  10.  Principal  component  scores  for  water  quality  PC1  at  five 
stations  on  the  Athabasca  River  during  1981. 


continued  . . 
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Figure  10.  Concluded. 
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Figure  11.  Principal  component  scores  for  water  quality  PC2  at  five 
stations  on  the  Athabasca  River  during  1981. 
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Figure  11.  Concluded 
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Figure  12.  Principal  component  scores  for  water  quality  PC3  at  five 
stations  on  the  Athabasca  River  during  1981. 
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Figure  13.  Principal  component  scores  for  water  quality  PC4  at  five 
stations  on  the  Athabasca  River  during  1981. 
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Figure  13.  Concluded 
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Figure  14.  Principal  component  scores  for  water  quality  PC5  at  five 
stations  on  the  Athabasca  River  during  1981. 
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Figure  14.  Concluded 
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Figure  15.  Principal  component  scores  for  water  quality  PC6  at  five 
stations  on  the  Athabasca  River  during  1981. 
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October.  Differences  between  left  and  right  banks  were  smaller  and 
less  consistent  at  stations  DA0207  and  DA0208.  As  was  the  case  with 
PCI , these  differences  between  the  left  and  right  sides  of  the 
Athabasca  River  appear  to  reflect  the  influences  of  right  bank 
tributary  streams,  particularly  the  Clearwater  River. 

PC3,  which  is  primarily  a measure  of  potassium,  lead,  zinc, 
and  copper,  generally  did  not  show  substantial  differences  between 
stations,  although  an  unusually  high  value  did  occur  near  the  first 
of  March  in  a left  bank  sample  at  Station  DA0207.  With  few  excep- 
tions, PC3  scores  were  similar  for  both  left  and  right  bank  samples 
at  all  stations. 

The  fourth  PC  (primarily  a measure  of  sodium,  chloride,  and 
ammonia)  had  much  higher  values  for  right  bank  samples  than  for  left 
bank  samples  at  Station  DA0203.  PC^  was  also  frequently  higher  for 
right  bank  samples  at  both  stations  DA0206  and  DA0207,  but  the  dif- 
ferences between  the  left  and  right  banks  were  not  as  large  or  as 
consistent  as  at  Station  DA0203.  These  differences  are  almost 
certainly  due  to  the  higher  sodium  and  chloride  concentrations 
present  in  right  bank  tributary  streams,  particularly  the  Clearwater 
River. 

PCS,  which  represents  primarily  dissolved  orthophosphate, 
showed  a general  declining  trend  from  high  values  in  the  spring  to 
low  values  in  the  fall,  with  some  fluctuations  in  between.  For  the 
most  part,  there  were  no  large  differences  between  stations,  but  PCS 
scores  were  frequently  higher  for  right  bank  samples  than  for  left 
bank  samples  at  stations  DA0203  and  DA0206. 

There  were  no  consistent  differences  between  left  and  right 
bank  samples  for  PC6,  which  is  primarily  a measure  of  phenolic  mater- 
ial. In  addition,  only  a few  differences  between  stations  were 
apparent.  Late  in  the  year  (after  September),  PC6  tended  to  be 
higher  at  stations  DA0206,  DA0207,  and  DA0208  than  at  Station 
DA0203.  Two  unusually  large  values  were  noted:  one  for  a mid-channel 
sample  taken  near  the  first  of  Dune  at  Station  CC0012,  and  one  for 
another  mid-channel  sample  taken  near  the  same  time  at  Station 
DA0206. 
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3.2  BENTHIC  INVERTEBRATES 

A separate  principal  component  analysis  was  performed  using 
benthic  invertebrate  abundance  data  from  each  of  the  months  of  May, 
3une,  3uly,  and  August  of  1981.  The  principal  components  so  derived 
were  then  used,  instead  of  the  individual  benthic  taxa,  for 
subsequent  analysis  of  variance  and  multiple  comparisons  tests 
comparing  sampling  stations.  Interpretation  of  principal  components 
was  based  simply  on  the  abundance  of  the  taxa  most  strongly 
associated  with  them.  For  reference  purposes,  the  abundance  data  for 
each  benthic  invertebrate  taxon  are  summarized  in  tables  20  to  23 
(Appendix  7.2). 

Although  the  benthic  invertebrate  data  were  subjected  to 
detailed  statistical  analysis,  the  reader  should  be  aware  that, 
because  of  the  way  in  which  benthic  sampling  was  conducted,  applica- 
tion of  statistical  procedures  is  not  entirely  valid.  Boerger  (1983) 
gives  the  following  description  of  the  sampling  procedure  used: 

At  each  station,  three  replicate  samples  were  collected, 
with  the  distance  between  samples  being  approximately 
10  m.  After  the  samples  were  placed  in  the  collecting  jar, 
they  were  examined  for  uniformity  of  sample  volume.  If  the 
volume  of  detritus  in  the  jars  differed  by  more  than  50%, 
more  samples  were  collected  and  the  three  most  similar  ones 
chosen.  This  procedure  was  based  on  the  correlation 
between  detritus  and  density  of  organisms  which  has  been 
frequently  reported  in  the  literature  (cf.  Rabeni  and 
Minshall  1977).  Such  a selection  procedure  should  reduce 
the  variability  between  samples  without  biasing  the  sample 
mean,  since  there  was  no  selection  for  either  high  or  low 
volumes  of  detritus,  but  only  for  a uniform  amount  of 
detritus. 

Given  that  the  density  of  benthic  invertebrates  is  cor- 
related with  detritus,  this  sample  selection  procedure  would  Indeed 
reduce  the  variability  among  samples,  the  result  being  that  sample 
variance  always  would  be  underestimated.  The  magnitude  of  the  bias 
is  unknown,  but  it  is  clear  that  it  might  be  large.  Since  it  appears 
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that  samples  were  collected  only  until  three  that  were  considered 
similar  enough  were  obtained,  the  sample  mean  cannot  be  expected  to 
be  the  same  as  the  mean  that  would  be  obtained  with  random  sampling 
(i.e.,  the  sample  mean  is  not  an  unbiased  estimate). 

The  effect  of  underestimated  variances,  in  terms  of  statis- 
tical hypothesis  testing,  would  be  a tendency  toward  accepting  dif- 
ferences as  being  statistically  significant  when,  in  fact,  they 
are  not.  Statistically  significant  differences  that  are  merely 
artifacts  of  the  sample  selection  procedure  would  tend  to  occur.  The 
effect  of  bias  of  the  sample  mean  is  less  predictable.  In  some 
instances  it  might  tend  to  create  false  significant  differences,  and 
in  others  it  might  tend  to  obscure  real  differences. 

In  view  of  the  above  factors,  which  seriously  affect  the 
validity  of  statistical  tests,  detailed  examination  of  all  signifi- 
cant differences  observed  among  stations  with  respect  to  benthic 
invertebrate  populations  was  considered  to  be  unwarranted.  The 
results  of  statistical  analyses  and  the  conclusions  drawn  from  them 
might  well  be  erroneous. 

In  the  following  sections,  the  results  of  analysis  of 
benthic  invertebrate  data  are  presented,  with  little  discussion,  for 
each  of  the  months  of  May,  June,  Ouly,  and  August.  This  is  followed 
by  a summary  discussion  that  is  restricted  to  consideration  of  only 
those  stations  that  might  be  expected  to  show  differences  related  to 
effluent  sources  or  to  the  influence  of  the  Clearwater  River. 

3.2.1  May  1981  Samples 

The  results  of  principal  component  analysis  of  benthic 
invertebrate  data  collected  in  May  are  presented  in  Table  3,  which 
contains  the  PC  loadings  matrix.  The  first  eight  principal  compon- 
ents accounted  for  61%  of  the  total  variance.  PCI  accounted  for 
1^.3%  of  the  total  variance,  and  the  amount  of  variance  explained  by 
each  of  PC3,  PC^,  PC5,  PC6,  PC7,  and  PCS  was  8.3%,  7.0%,  5.9%,  7.2%, 
6.6%,  6.2%,  and  5.6%,  respectively.  Details  of  which  benthic 
invertebrate  taxa  are  most  closely  associated  (i.e.,  have  loadings  of 
at  least  0.25  absolute  value)  with  each  PC  are  outlined  in  Table 
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Principal  component  scores  for  each  benthic  sample  are  determined 
primarily  by  the  abundances  of  the  taxa  indicated  in  Table  4 as  being 
associated  with  each  PC. 

Analysis  of  variance,  which  was  performed  for  each  PC 
(Table  5),  indicated  that  there  were  significant  differences  among 
sampling  stations  only  with  respect  to  PC2,  PC3,  PC6,  and  PC7. 
Multiple  comparisons  tests  (which  make  all  pairwise  comparisons 
between  stations)  for  each  of  these  PCs  are  summarized  in  Table  6. 
The  tests  are  for  differences  between  mean  PC  scores,  but  the  mean 
abundances  for  the  benthic  taxa  most  closely  associated  with  each  PC 
have  been  included  to  facilitate  relating  differences  in  PC  scores  to 
the  specific  taxa  involved  in  each  case.  There  were  no  significant 
differences  in  PC7  scores  among  station  pairs  despite  the  fact  that 
analysis  of  variance  detected  significant  (P  < 0.05)  variation  in  PC7 
scores.  This  apparent  contradiction  is  a result  of  the  lower  sensi- 
tivity of  the  multiple  comparisons  procedure. 

3.2.2  June  1981  Samples 

The  principal  component  loadings  matrix  for  benthic  inver- 
tebrate samples  collected  during  the  month  of  June  is  presented  in 
Table  7.  In  this  sampling  period,  seven  principal  components 
accounted  for  58%  of  the  total  variance  in  benthic  invertebrate  abun- 
dances. PC1  accounted  for  16.4%  and  PC2  accounted  for  11.6%  of  the 
total  variance,  while  the  amount  of  variance  explained  by  each  of  PCs 
3,  4,  5,  6,  and  7 was  7.0%,  7.4%,  6.2%,  5.3%,  and  4.3%,  respec- 
tively. The  interpretations  of  principal  components,  in  terms  of  the 
benthic  taxa  most  closely  associated  with  each  PC,  are  described  in 
Table  8. 

Analysis  of  variance  results  comparing  sampling  stations 
with  respect  to  PC  scores  are  summarized  in  Table  9.  Significant 
differences  among  stations  occurred  only  for  PC2,  PC3,  and  PC4.  The 
specific  differences  between  stations,  with  respect  to  these  three 
PCs,  were  determined  using  the  Student-Newman-Keuls  multiple  compari- 
sons procedure,  and  the  results  of  these  comparisons  are  summarized 
in  Table  10. 
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Table  4.  Benthic  invertebrate  taxa  with  the  highest  absolute  PC 

loadings  (at  least  0.25),  May  samples.  Principal  component 
scores  are  determined  primarily  by  the  log  (no./m^  + 1) 
of  these  taxa. 


PC1: 

Heptagenia 

0.663 

Metretopus 

-0.624 

Perlodidae 

0.633 

vs. 

Siphlonurus 

-0.487 

Rhithrogena 

0.438 

Oligochaeta 

-0.291 

Chironomidae 

0.379 

Simuliidae 

-0.265 

Ephemerella 

0.324 

Baetisca 

0.279 

Empididae 

0.271 

Baetis  sp.  A 

0.258 

PC2: 

Oligochaeta 

0.833 

Chironomidae 

0.556 

Empididae 

0.373 

Nematomorpha 

0.321 

PC3: 

Chironomidae 

0.633 

vs. 

Oligochaeta 

-0.395 

Perlodidae 

-0.389 

Simuliidae 

-0.369 

Heptagenia 

-0.330 

Hydropsychidae 

-0.255 

PC4: 

Metretopus 

0.658 

vs. 

Empididae 

-0.355 

Heptagenia 

0.458 

PC5: 

Simuliidae 

0.692 

vs. 

Ephemerella 

-0.318 

Corixidae 

0.379 

Ceratopogonidae 

-0.277 

Chironomidae 

0.286 

Empididae 

-0.259 

continued  . . . 
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Table  4.  Concluded. 

PC6:  Empididae  0.515 

Perlodidae  0.456 

Simuliidae  0.434 

Siphlonurus  0.301 

Nematomorpha  0.280 


PC7; 

Ephemerella 

0.602 

vs. 

Ophiogomphus 

-0.464 

Siphlonurus 

-0.426 

PCS: 

Nematomorpha 

0.577 

vs. 

Baetis  sp.  A 

-0.401 

Empididae 

0.267 

Hydropsychidae 

-0.315 

Heptagenia 

0.264 

Perlodidae 

-0.290 
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Table  5.  Analysis  of  variance  of  benthic  invertebrate  principal 

component  scores  at  the  16  sampling  stations,  May  samples. 


Between  Stations 

Error 

Principal 

Component 

Degrees 

of  Sum  of 

Freedom  Squares 

Degrees 
Mean  of 

Squares  Freedom 

Sum  of  Mean  F 

Squares  Squares 

PCI 

15 

12.686 

0.846 

80 

82.314 

1.029 

0.82 

PC2 

15 

34.152 

2.277 

80 

60.848 

0.761 

2,99*** 

PC3 

15 

32.118 

2.141 

80 

62.882 

0.786 

2.72** 

PC4 

15 

21.763 

1.451 

80 

73.237 

0.916 

1.59 

PC5 

15 

23.625 

1.575 

80 

71.376 

0.892 

1.77 

PC6 

15 

36.688 

2.446 

80 

58.312 

0.729 

3.36*** 

PC7 

15 

24.612 

1.641 

80 

70.387 

0.880 

1.87* 

PCS 

15 

10.925 

0.728 

80 

84.075 

1.051 

0.69 

* significant  with  P < 0.05 

**  significant  with  P < 0.01 

***  significant  with  P < 0.001 
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Table  7.  Principal  component  loadings  for  the  first  seven  principal  components  of  benthic  invertebrate 
abundance,  Oune  samples.  Loadings  greater  than  0.25  absolute  value  are  printed  in  boldface. 
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Table  8 


PC1: 


PC2: 


Benthic  invertebrate  taxa  with  the  highest  absolute  PC 
loadings  (at  least  0.25),  Gune  samples.  Principal 
component  scores  are  determined  primarily  by  the  log 
(no./m^  + 1)  of  these  taxa. 


Chironomidae 

0.747 

vs.  Perlodidae 

-0.291 

Tricorythodes 

0.738 

Ephemerella 

-0.290 

Caenis 

0.719 

Isonychia 

6.671 

Heptagenia 

0.614 

Cloeon  sp.  1 

0.503 

Oligochaeta 

0.438 

Ophiogomphus 

0.421 

Hydropsychidae 

0.370 

Brachycercus 

0.366 

Simuliidae 

0.327 

Pseudocloeon 

0.278 

Perlodidae 

0.749 

vs.  Oligochaeta 

-0.476 

Baetis  sp.  A 

0.696 

Simuliidae 

-0.266 

Heptagenia 

0.559 

Ametropus 

-0.250 

Ephemerella 

0.484 

Metretopus 

0.399 

Pseudocloeon 

0.393 

Hydracarina 

0.276 

Ophiogomphus 

0.271 

continued  . . 
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Table  8.  Concluded. 


PC3; 

Oligochaeta 

0.584 

vs.  Ameletus 

-0.378 

Baetis  sp.  A 

0.422 

Cloeon  sp.  1 

-0.377 

Simuliidae 

0.416 

Caenis 

-0.373 

Nematomorpha 

0.342 

Ceratopogonidae 

0.303 

Isonychia 

0.253 

PC4: 

Ceratopogonidae 

0.481 

vs.  Hydropsychidae 

-0.500 

Caenis 

0.375 

Ameletus 

-0.406 

Brachycercus 

0.362 

Heptaqenia 

-0.346 

Ophioqomphus 

0.337 

Tricorythodes 

-0.298 

Pseudocloeon 

0.316 

Cloeon  sp.  1 

0.299 

PC5; 

Ametropus 

0.677 

vs.  Brachycercus 

-0.348 

Pseudocloeon 

0.373 

Baetis  sp.  A 

-0.281 

Ephemerella 

0.351 

Tricorythodes 

0.309 

Centroptilum  sp.  1 

0.305 

Ceratopogonidae 

0.299 

PC6: 

Metretopus 

0.620 

vs.  Ceratopogonidae 

-0.543 

Simuliidae 

-0.334 

PC7; 

Chironomidae 

0.359 

vs.  Metretopus 

-0.399 

Ephemerella 

0.310 

Isonychia 

-0.292 

Pseudocloeon 

0.307 

Baetis  sp.  C 

0.272 

Ophioqomphus 

0.267 
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Table  9.  Analysis  of  variance  of  benthic  invertebrate  principal 

component  scores  at  the  16  sampling  stations,  June  samples. 


Between  Stations  Error 


Principal 

Component 

Degrees 

of  Sum  of 

Freedom  Squares 

Mean 

Squares 

Degrees 

of 

Freedom 

Sum  of  Mean 
Squares  Squares 

F 

PC1 

15 

22.^72 

1.498 

80 

72.529 

0.907 

1.65 

PC2 

15 

55.091 

3.673 

80 

39.909 

0.499 

7.36*** 

PC3 

15 

35.367 

2.358 

80 

59.634 

0.745 

3.16*** 

PC4 

15 

^8.^86 

3.232 

80 

46.514 

0.581 

5.56*** 

PC5 

15 

23.612 

1.574 

80 

71.389 

0.892 

1.76 

PC6 

15 

21.217 

1.414 

80 

73.783 

0.922 

1.53 

PC7 

15 

14.673 

0.978 

80 

80.327 

1.004 

0.97 

***  significant 

with  P < 0. 

.001 
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Table  10.  Concluded 
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3.2.3  July  1981  Samples 

The  principal  component  loadings  for  benthic  taxa  collected 
during  the  month  of  3uly  are  presented  in  Table  11.  In  Ouly,  the 
first  eight  principal  components  accounted  for  65%  of  the  total  vari- 
ance in  the  abundances  of  the  taxa  collected.  The  amount  of  variance 
explained  by  each  of  PCI,  PC2,  and  PC3  was  15.5%,  13.3%,  and  8.9%. 
PC^,  PC5,  PC6,  PC7,  and  PCS  accounted  for  6.0%,  6.7%,  5.A-%,  5.3%,  and 
4.3%,  respectively,  of  the  total  variance.  The  composition  of  the 
eight  principal  components,  in  terms  of  the  major  constituent  taxa, 
are  described  for  each  PC  in  Table  12. 

Analysis  of  variance  of  the  principal  component  scores  for 
each  PC  indicated  that  significant  differences  among  stations 
involved  all  principal  components  except  PCS.  Multiple  comparisons 
tests,  which  compare  all  pairs  of  stations,  were  performed  for  all 
PCs  for  which  a significant  difference  was  indicated  by  analysis  of 
variance.  The  results  of  these  pairwise  comparisons  are  summarized 
in  Table  14.  In  3uly,  there  were  more  significant  differences 
between  stations,  involving  more  principal  components,  than  in  either 
May  or  June. 

3.2.4  August  1981  Samples 

Results  of  principal  component  analysis  of  the  August  ben- 
thic invertebrate  samples  are  presented  in  Table  15,  which  includes 
the  PC  loadings  for  each  of  the  benthic  taxa.  The  first  eight  prin- 
cipal components  accounted  for  68%  of  the  total  variance  of  the  ben- 
thic invertebrate  abundances  in  August.  PCI  accounted  for  15.0%  and 
PC2  accounted  for  12.9%  of  the  total  variance.  PCs  3,  4,  5,  6,  7, 
and  8 explained  8.9%,  7.2%,  7.3%,  5.8%,  5.7%,  and  5.0%,  respectively, 
of  the  total  variance.  Descriptions  of  the  eight  principal  compon- 
ents, in  terms  of  the  benthic  taxa  most  closely  associated  with  each 
PC,  are  given  in  Table  16.  Five  of  the  principal  components  (PC1 , 
PC2,  PC3,  PC4,  and  PC6)  showed  significant  differences  among  stations 
when  tested  by  analysis  of  variance  (Table  17).  The  details  of  which 
stations  were  significantly  different  from  which  other  stations  (as 
determined  by  multiple  comparisons  tests)  are  summarized  in  Table  18. 


Table  11.  Principal  component  loadings  for  the  first  eight  principal  components  of  benthic  invertebrate 
abundance,  3uly  samples.  Loadings  greater  than  0.25  absolute  value  are  printed  in  boldface. 
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Table  1 


PC1: 


PC2: 


PC3; 


. Benthic  invertebrate  taxa  with  the  highest  absolute  PC 
loadings  (at  least  0.25),  3uly  samples.  Principal 
component  scores  are  determined  primarily  by  the  log 
(no./m^  + 1)  of  these  taxa. 


Caenis 

0.829 

Tricorythodes 

0.7U 

Heptagenia 

0.677 

Ophiogomphus 

0.588 

Isonychia 

0.^75 

Cladocera 

0.439 

Cloeon  sp.  1 

0.429 

Collembola 

0.382 

Centroptilum  sp.  1 

0.334 

Ephemerella 

0.327 

Ameletus 

0.301 

Nematomorpha 

0.277 

Cloeon  sp.  1 

0.757 

vs.  Isonychia 

-0.535 

Cladocera 

0.734 

Hydropsychidae 

-0.533 

Chironomidae 

-0.490 

Tricorythodes 

-0.436 

Baetis  sp.  X 

-0.388 

Heptagenia 

-0.383 

Ephemerella 

-0.310 

Oligochaeta 

-0.281 

Nematomorpha 

0.769 

vs.  Pseudocloeon 

-0.473 

Oligochaeta 

0.691 

Heptagenia 

-0.301 

Baetis  sp.  A 

0.402 

Baetis  sp.  X 

-0.264 

Ostracoda 

0.279 

Chironomidae 

0.269 

continued  . . . . 
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Table  12.  Concluded. 


PC4: 

Ametropus 

0.598 

vs.  Hydropsychidae 

-0.359 

Ceratopogonidae 

0.317 

Isonychia 

-0.351 

Ameletus 

-0.326 

Ephemerella 

-0.315 

Cloeon  sp.  1 

-0.253 

PC5: 

Ostracoda 

0.506 

vs.  Baetis  sp.  A 

-0.409 

Chironomidae 

0.478 

Isonychia 

-0.369 

Collembola 

0.388 

Ametropus 

-0.281 

Oligochaeta 

0.337 

Pseudocloeon 

0.312 

Leptophlebia 

0.311 

PC6: 

Collembola 

0.478 

vs.  Hydropsychidae 

-0.332 

Pseudocloeon 

0.411 

Ephemerella 

-0.322 

Ophiogomphus 

0.401 

Ostracoda 

-0.316 

Baetis  sp.  A 

0.261 

PC7; 

Ameletus 

0.619 

vs.  Baetis  sp.  A 

-0.428 

Ametropus 

0.449 

Collembola 

-0.274 

Ophiogomphus 

0.304 

PCS: 

Baetis  sp.  X 

0.707 

vs.  Ametropus 

-0.369 

Oligochaeta 

0.266 

Hydropsychidae 

-0.299 
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Table  13.  Analysis  of  variance  of  benthic  invertebrate  principal 

component  scores  at  the  16  sampling  stations,  Duly  samples. 

Between  Stations  Error 


Degrees  Degrees 

Principal  of  Sum  of  Mean  of  Sum  of  Mean  F 

Component  Freedom  Squares  Squares  Freedom  Squares  Squares 


PCI 

15 

53.456 

3.564 

80 

41.544 

0.519 

6.86*** 

PC2 

15 

31 .497 

2.100 

80 

63.504 

0.794 

2.65** 

PC3 

15 

33.068 

2.205 

80 

61.933 

0.774 

2.85** 

PC4 

15 

39.868 

2.658 

80 

55.131 

0.689 

3.86*** 

PC5 

15 

38.159 

2.544 

80 

56.841 

0.711 

3.58*** 

PC6 

15 

28.610 

1.907 

80 

66.391 

0.830 

2.30** 

PC7 

15 

27.995 

1.866 

80 

67.005 

0.838 

2.23* 

PCS 

15 

20.405 

1.360 

80 

74.595 

0.932 

1.46 

* significant 

with 

P 

< 

0.05 

**  significant 

with 

P 

< 

0.01 

***  significant 

with 

P 

< 

0.001 
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Table  15.  Principal  component  loadings  for  the  first  eight  principal  components  of  benthic 

invertebrate  abundance,  August  samples.  Loadings  greater  than  0.25  absolute  value  are 
printed  in  boldface. 
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Table 


PCI: 


PC2: 


16.  Benthic  invertebrate  taxa  with  the  highest  absolute  PC 
loadings  (at  least  0.25),  August  samples.  Principal 
component  scores  are  determined  primarily  by  the  log 
(no./m^  + 1)  of  these  taxa. 


Heptagenia 

0.622 

Ceratopogonidae 

0.563 

Tricorythodes 

0.560 

Cloeon  sp.  1 

0.537 

Baetis  sp.  X 

0.528 

Empididae 

0.^99 

Oligochaeta 

0.487 

Hydropsychidae 

0.458 

Cladocera 

0.431 

Ostracoda 

0.418 

Perlodidae 

0.396 

Ametropus 

0.368 

Ophiogomphus 

0.367 

Neureclipsis 

0.285 

Caenis 

0.279 

Perlodidae 

0.700 

vs.  Chironomidae 

-0.527 

Hydropsychidae 

0.678 

Leptophlebia 

-0.524 

Ametropus 

0.547 

Neureclipsis 

-0.472 

Caenis 

0.312 

Ostracoda 

-0.472 

Metretopus 

0.283 

Nematomorpha 

-0.350 

Cladocera 

-0.320 

Cloeon  sp.  1 

-0.306 

continued  .... 
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Table  16. 

Continued. 

PC3: 

Oligochaeta 

0.767 

vs.  Leptophlebia 

-0.463 

Chironomidae 

0.443 

Ophioqomphus 

-0.445 

Nematomorpha 

0.369 

Heptaqenia 

-0.370 

Empididae 

0.336 

Baetis  sp.  X 

-0.369 

Cladocera 

-0.326 

PC^: 

Baetis  sp.  X 

0.539 

vs.  Empididae 

-0.486 

Nematomorpha 

0.365 

Cloen  sp.  1 

-0.450 

Pseudocloeon 

0.358 

Metretopus 

-0.361 

Perlodidae 

0.353 

Leptophlebia 

0.302 

PC5: 

Baetis  sp.  C 

0.505 

vs.  Ceratopogonidae 

-0.556 

Pseudocloeon 

0.470 

Metretopus 

-0.331 

Ophioqromphus 

0.333 

Ostracoda 

-0.260 

Nematomorpha 

0.324 

Tricorythodes 

0.314 

Hydropsychidae 

0.304 

Neureclipsis 

0.263 

PC6: 

Empididae 

0.384 

Tricorythodes 

-0.413 

Oecetis 

0.342 

Ametropus 

-0.375 

Perlodidae 

0.337 

Baetis  sp.  X 

-0.310 

Neureclipsis 

0.305 

Caenis 

-0.310 

Ceratopogonidae 

0.275 

Heptaqenia 

0.268 

continued  .... 
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Table  16. 

Concluded. 

PC7: 

Neureclipsis 

0.394 

vs.  Chironomidae 

-0.346 

Centroptilum  sp.  1 

0.380 

Hydropsychidae 

-0.303 

Ophiogomphus 

0.352 

Cloeon  sp.  1 

-0.251 

Ametropus 

0.326 

Cladocera 

0.315 

Oecetis 

0.269 

Caenis 

0.260 

PCS: 

Cladocera 

0.555 

vs.  Tricorythodes 

-0.341 

Metretopus 

-0.339 

Baetis  sp.  C 

-0.319 

Pseudocloeon 

-0.285 

Ceratopogonidae 

-0.274 
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Table  17.  Analysis  of  variance  of  benthic  invertebrate  principal 
component  scores  at  the  16  sampling  stations,  August 
samples. 


Between  Stations  Error 


Principal 

Component 

Degrees 

of  Sum  of 

Freedom  Squares 

Mean 

Squares 

Degrees 

of 

Freedom 

Sum  of  Mean 
Squares  Squares 

F 

PCI 

15 

29.051 

1.937 

32 

17.949 

0.561 

3.45** 

PC2 

15 

35.924 

2.395 

32 

11.076 

0.346 

5,92*** 

PC3 

15 

29.803 

1.987 

32 

17.197 

0.537 

3.70*** 

PC4 

15 

29.479 

1.972 

32 

17.421 

0.544 

3.62** 

PC5 

15 

20.196 

1.346 

32 

26.804 

0.838 

1.61 

PC6 

15 

31.769 

2.118 

32 

15.231 

0.476 

Zf , 45*** 

PC7 

15 

19.903 

1.327 

32 

27.098 

0.847 

1.57 

PCS 

15 

18.408 

1.227 

32 

28.592 

0.894 

1.37 

**  significant  with  P < 0.01 

***  significant  with  P < 0.001 


60 


4->  to 
C 4) 

(0  JZ  JZ 
0 0)0 
•H  .H  fO 
<t-  JZ  O 


-D  ^ 


C O tt- 
0>JZ  O 


> *H 
C 4_ 


4-> 

fO 

Q. 


•H 

•H  JZ 

ft) 

0 

C 4->  4-> 

X 

Q. 

o 

O)  fO 

ft) 

Ql 

•H 

•H  >»JZ 

JJ 

Z) 

^ CO  ^ 4J 

4-> 

t<- 

0 

c 

-D  to 

o 

SI 

o 

• O C 

XJ 

to  c o 

to 

O *H  *H 

4) 

c 

t>- 

'-H  O 

O 

•H 

o 

Q.  -h-} 

c 

E 4-> 

ft) 

*D 

to 

fO  <0  to 

*D 

0 

to  c 

C 

•o 

to 

o >» 

O 

D 

4) 

4J  *H  I-H 

-Q 

I-H 

4J 

to  4-»  C 

< 

0 

O ft)  o 

C 

to 

O)  4J 

•H 

c 

O CO 

• 

o 

< 

•o 

0 

to 

u 

4) 

•H 

^ . 4) 

•o 

ft) 

(-1 

C to  JZ 

o 

o O +J 

1— 1 

O 

Q.*H  O 

O Q_ 

E 

^ o 

c 

O 

ft)  O jZ 

•H 

JZ 

o 

JJ  to  o 

o 

to  fO 

4) 

ft) 

4) 

O 4) 

Cj 

0 

f— 1 

0)0. 

ft) 

Q. 

C E 

c 

•H 

♦H  C O 

C 

o 

^ fO  fcj 

o 

i-H 

CL  4)  t(_ 

iH 

2 

E E 

XJ 

0 

E 

fO 

ft) 

D 

to  4)  LA 

x-> 

I-H 

to 

JZ  . 

to 

ft) 

f— 1 

JZ  +->  o 

> 

ZJ 

O 

4; 

ft)  ^ 

4) 

0 

4)  O 

JZ 

X-> 

I-H  Q_ 

XJ 

D 

1 

JJ)  4) 

o 

I-H 

ft)  -Q 

o 

c 

4J» 

4) 

to 

ft) 

to  to  C 

C 

E 

4)  4)  4) 

o 

ft) 

? 

t-i  C C-i 

4> 

O *H  4) 

XJ 

LA 

Z 

O ^ <f- 

to 

AJ 

tO  tt_ 

ft) 

1 

4)  *H 

4)  O 

JH  JZ  "D 

I-H 

4-> 

C 4J 

J-> 

C 

4)  >,  JH 

to 

o 

C 

fO 

ft) 

•O 

O J2  j-> 

0 

D 

Q.  C 

E 

I-H 

+-> 

E T3  ft) 

O 

CO 

O 4)  O 

(-1 

X-» 

O 4J  *H 

<4- 

ft) 

t»- 

ft)  4- 

o 

i-H  O *H 

XJ 

ft)  -H  C 

c 

to 

to 

Q.-0  O) 

4) 

O) 

4-> 

•H  C *H 

Sx 

c 

I-H 

O *H  to 

0 

•H 

o 

c 

<L- 

■D  ' 

<0 

•H  4)  HJ 

<4- 

ft)  . 

o 

t-i  u.  o 

•H 

O 

cc 

CL  f^  C 

■D 

I-H 

O fn  O rri 


VO  O ro  ( 

m o Jl-  ' 

^ fvi 


O O O 

o VO  \o  d o 


r*.  o o o 


o o o r*  o 


mOf^OvoOO'O'OO 
(M  VO  VO  (M 


f^OOr^rrvr^rrvoOOr^'^f^OO 
vodd 


o o o 

d d d 


oovor'»oo<^o<^oo 

d d 


o o o o 


o o f'-  r"  o 
O O VO  rrv  VO  VO  O 


VO 

n^mOOf^OOOOOOOOOO  vo 

rrvrriOO'^OOOOOOOOOO  <- 
VO  ' 

ao 


fvj(»vOf^ooovr^oov— oovocor^inCv 
cjvo'O'^'MCTvcoiAfnv— ovvovooor^ 

voiAvAiAi/v^^<j--^^(^fnr>^(Mfvj 

ddddddodddddooo 


c 

a; 

crt  o 


Q.  a.  V <; 


« (0  (O  (O 


^ 0) 
o (0  ifl  (T3 

>v  t-  -o  -D 

in  aj  o — 
Q.  o o "a 


o.  u 

V V 

X o 


•o  j->  ti  aj 


O X O O Q.  < 


continued 


Mean  No./m^  at  Each  Station 


81 


■o 

<D 

D 

C 


o o o o o i 
d o o o o 


r^or'•^^omor^OO(»^ 
d VO  d d d 


r^r^OOOOOOOOOf^ 

dvodddddddddd 


OrrvfnoOOOOOr'-Or' 

dddddddddvoovo 

O O lA  fVJ  (M 

M o 

<nr*.(T\oor>.Of^OO'Ar^ 

mvomoovoofAooiAvo 


r^fAmor^fAor'Oooo 

vOfAiAOvOrrvOvOOOOO 


Or^OOOOOOOfAOO 

ovoddddddoddo 


r'0'Aoooooomor» 

'OOi^OOOOOO'AO'O 
«-  VO 

o 

Or^r^ooofAoor^r^O 


OiAP^mrrir^f^OO 


o o o o o 
o d VO  d o d 

VO 

00  Ov 


lA 

I O VO 

i d d 


rrvOfAf^Or^P^fAO 

moiAvoovovorAO 

»-  T- 

(A  O V- 


o rv 
d VO 


OOiAOfACOO 

OOiAOfAOOO 
OV  AJ 

o 


r'Or^OO'AiArAO 

vOOvOOOiArrvmo 
fvj  lA  (M  rvj  ^ Ov  <A 
(AO  <M 


O O O O O O O O O !>.  O O 

ddddoddddvodd 


OOt^(A(AOOOA' 

OOvOfAfAOOOvO 


r*.  o (A  ( 

VO  O (A  I 


r'  O rA 


O O O O 


t>.ooO(ArAt>.r^or^o 


o o o o 


0(AOOO(A(AvovorAOiA 
VA  VO  (M  VO 
Ov  V- 


fA  o fA 


O (A 
«M  lA 

AJ  Ov 


OflOr^AJ(A»>.^fVJAJOOVO 

Or'«^r-®(MAjr^l>.lAAJO 

r^vOIArAAJl/VlA^^(A(A(A 

dddddddddddd 


t^(AOvvO(AiAOOvvo 

VO^VO«A\p^r^VOAJ 

f'.<f(A(A^^rArA(A 

OOOOCOOOO 


■O  ^ 

•g  ^ 
■ O JZ 


j: 

Q. 

ki  (O  Q. 

■a  o t-  <fl 

O E lU 

o o u c 

^ 4J  O O 
(O  TO  D 

e T)  o 

(O  4)  -H 

* ^ u o 


1 s 


a.  ^ 
ki  1)  o> 
O (0  -H 
E -o  j: 


E Q. 

4)  E 4) 
Z Ld 


continued 


Table  18.  Concluded 


82 


■ o o o o o o 

• rn  O O O O O O 


OvOrnof^r^OO 


O O 

O d 


O O O 

d d d d 


O 

O VO  va  VO  VO 


'^oor^or^or^ 


ooor'Orvor^ 

OOOvOOvOOve 


ovirvoo(n«Mvoo»- 

(^voi/virvOccvrvvo 

mrrvf»vf»v(TVtft^fn 

ododdddd 


• </} 
O.  3 
10  a. 


u a.  CL  o 


I 

I 


83 


3.2.5  Summary  of  Benthic  Invertebrate  Analyses 

Is  is  clear  from  Tables  6,  10,  14,  and  18  that  most 
stations  differed  significantly  from  others  at  lease  once  with 
respect  to  at  least  one  principal  component.  In  the  discussion  to 
follow,  however,  only  those  differences  are  discussed  that  may  be 
attributable  to  (1)  Fort  McMurray  sewage  discharge,  (2)  Clearwater 
River  influences  and  (3)  Suncor  plant  effluent  (i.e..  Station  7 is 
compared  with  Station  8,  and  Station  4 is  compared  with  Station  5). 
It  must  be  remembered  throughout  that  differences  also  may  be  simply 
artifacts  of  the  sample  selection  method  used. 

3. 2. 5.1  Stations  7E,  7W,  8E,  and  8W.  Stations  8E  and  8W  are 
controls  above  the  influences  of  Fort  McMurray  and  the  Clearwater 
River.  Station  7E  is  within  the  plume  of  the  Clearwater  River,  and 
Station  7W  is  thought  to  be  influenced  by  the  Fort  McMurray  sewage 
treatment  plant  effluent  (Boerger  1983). 

No  differences  in  the  benthic  invertebrate  principal 
components  for  these  stations  were  detectable  in  May  (Table  6);  how- 
ever several  significant  differences  were  detected  among  the  stations 
in  the  Dune  data  (PC2  and  PC4;  Table  10).  Station  7W  had  a signifi- 
cantly higher  Dune  PC2  score  than  8E,  primarily  because  of 
higher  numbers  of  Heptaqenia  at  7W.  Station  7W  did  not  differ 
significantly  from  Station  8W  with  respect  to  Dune  PC2,  but  Station 
7E  had  a significantly  higher  PC2  score  than  both  stations  8E  and 
8W.  The  difference  was  largely  due  to  greater  numbers  of  certain 
Ephemeroptera  (Heptaqenia,  Baetis  sp.  A,  Metretopus,  Pseudocloeon) 
and  fewer  Oligachaeta  at  Station  7E.  Station  7E  did  have  a signifi- 
cantly higher  Dune  PC4  score  than  Station  7W,  primarily  because  of 
the  much  greater  abundance  of  Caenis  at  Station  7E. 

In  Duly,  differences  among  the  four  stations  were  confined 
to  PCI  and  PC5  (Table  14).  Stations  7E  and  7W  had  higher  PCI  values 
than  stations  8E  and  8W,  principally  because  certain  Ephemeroptera 
(especially  Caenis,  Tricorythodes,  and  Heptaqenia)  and  (in  the  case 
of  Station  7E)  Cladocera  and  Cloeon  were  more  abundant  at  the  two 


downstream  stations.  Station  7W  had  a much  higher  PC5  score  than  the 
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other  three  stations  because  of  its  abundant  population  of 
Oligochaeta. 

More  differences  among  the  benthic  invertebrate  associa- 
tions of  the  four  stations  were  found  in  August,  with  respect  to  PC1 , 
PC2,  and  PC3  scores  (Table  18).  Station  7E  had  a significantly 
higher  PC1  score  than  Station  8E,  primarily  because  of  the  greater 
numbers  of  certain  Ephemeroptera  (especially  Heptagenia, 
Tricorythodes , Cloeon  sp.  1,  and  Baetis  sp.  X),  Ceratopogonidae, 
Cladocera,  Ostracoda,  Ophiogomphus,  and  Neureclipsis.  Station  7E 
also  had  a lower  August  PC2  score  than  stations  8W  and  7W.  The 
difference  was  caused  mainly  by  the  low  numbers  of  Perlodidae  and 
Hydropsychidae  and  the  greater  numbers  of  Oligochaeta,  Leptophlebia, 
Neureclipsis,  Ostracoda,  Cladocera,  and  Cloeon  sp.  1 at  Station  7E. 
Finally,  stations  7E  and  7W  both  had  lower  August  PC3  scores  than 
stations  8E  and  8W,  mostly  due  to  various  combinations  of  lower 
numbers  of  Oligochaeta  and/or  Chironomidae  and  higher  numbers  of 
various  positive-loading  taxa  at  Stations  7E  and  7W. 

It  is  interesting  to  note  that  stations  8E  and  8W  did  not 
differ  significantly  from  each  other  with  respect  to  the  31  principal 
components  of  the  benthic  invertebrate  abundance  data  over  the  entire 
study  period.  It  would  appear  that  these  stations  are  sufficiently 
alike  in  their  benthic  invertebrate  associations  that  either  or  both 
could  serve  as  control  sites. 

Station  7W,  thought  to  be  influenced  by  Fort  McMurray 
sewage  discharge  (Boerger  1983),  might  be  expected  to  show  evidence 
of  enrichment  such  as  increased  abundances  of  algae-eating  inverte- 
brates. This  in  fact  did  account  for  differences  between  certain 
principal  component  scores  at  Station  7W  and  the  control  stations. 
Most  of  the  Ephemeroptera  (except  Metretopus,  thought  to  be  preda- 
ceous by  Merritt  and  Cummins  1978),  whose  greater  abundance  at 
Station  7W  accounted  for  certain  significantly  higher  principal 
component  scores  there,  are  at  least  partly  algae-scrapers  (Merritt 
and  Cummins  1978).  However,  some  evidence  is  inconsistent.  Station 
7W  had  a significantly  increased  July  PCS  score  primarily  because  of 
high  numbers  of  Oligochaeta  there.  The  situation  was  reversed  in 
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August,  when  very  low  numbers  of  Oligochaeta  at  Station  7W  gave  it  a 
PC3  score  significantly  lower  than  that  at  one  of  the  control 
stations. 

Station  7E,  in  the  Clearwater  River  plume,  differed  from 
the  control  stations  with  regard  to  several  of  the  principal  compon- 
ents. Some,  but  not  all,  of  the  differences  were  due  principally  to 
higher  numbers  of  certain  taxa  at  Station  7E  that  are  characteristic 
of  slow-flowing  or  even  lentic  environments,  such  as  Caenis,  Metre- 
topus,  Cladocera,  Cloeon,  Ceratopogonidae,  Neureclipsis,  and 
Leptophlebia  (Merritt  and  Cummins  1978;  Pennak  1978).  Water  velocity 
at  Station  7E  was  not  clearly  slower  than  at  least  one  of  the  control 
stations  (Boerger  1983).  It  is  possible  that  numbers  of  these 
characteristically  lentic  or  slow-water  groups  were  contributed  by 
drift  out  of  the  Clearwater  River. 

3. 2. 5. 2 Stations  4E,  ^W,  3E,  and  3W.  Stations  5E  and  5W  are 
located  a short  distance  upstream  of  the  Suncor  development,  and  were 
evidently  intended  to  act  as  control  sites.  Unfortunately,  as 
Boerger  (1983)  has  pointed  out,  there  are  no  suitable  control  sites 
upstream  of  Suncor  because  waters  from  the  Clearwater  River  and  the 
Athabasca  River  flow  separately  and  unmixed  past  5E  and  5W  until  the 
river  is  forced  into  a sharp  turn  by  the  embankment  at  the  upstream 
end  of  the  Suncor  site.  Presumably,  benthic  invertebrate  associa- 
tions at  Station  5E  reflect  the  influence  of  Clearwater  water,  and 
those  at  Station  5W  reflect  conditions  in  the  Athabasca  water.  In 
fact,  the  two  control  stations  differed  significantly  with  respect  to 
three  principal  components  of  the  benthic  invertebrate  data  (May  PC3, 
June  PC4,  and  3uly  PCI;  Tables  6,  10  and  1^,  respectively).  Inverte- 
brate associations  at  stations  4E  and  4W,  on  the  other  hand,  reflect 
the  influence  of  the  mixed  Clearwater  and  Athabasca  water  as  well  as 
the  influence  of  the  Suncor  plant,  if  any.  In  addition.  Station  5W 
is  located  downstream  from  Poplar  Creek,  through  which  saline  mine 
depressurization  water  from  Syncrude  Canada  Ltd.  is  discharged. 

None  of  the  31  principal  components  examined  showed 
significant  differences  between  either  of  the  Station  4 sites  and 
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both  the  control  sites  (Tables  6,  10,  1^,  and  18).  To  this  extent, 
there  is  no  evidence  that  Suncor  operations  have  had  a catastrophic 
effect  on  the  river  macrobenthos. 

3.3  RELATIONSHIPS  BETWEEN  WATER  QUALITY  AND  BENTHIC 

INVERTEBRATES 

In  1981,  water  quality  data  were  collected  for  only  five  of 
the  eight  benthic  sampling  stations  (stations  1,  2,  3,  6,  and  8). 
Examination  of  correlations  between  water  quality  and  benthic  inver- 
tebrates is  therefore  limited  to  these  five  stations.  The  results  of 
correlation  analyses  of  benthic  invertebrate  principal  component 
scores  and  water  quality  principal  component  scores  are  presented  in 
Table  19.  Fourteen  correlation  coefficients  were  found  to  be 
significant  at  P < 0.05,  11  were  significant  at  P < 0.01,  and  12  were 
significant  at  P < 0.001. 

In  the  discussion  to  follow,  only  correlations  significant 
at  P < 0.001  are  considered.  The  large  number  of  individual  analyses 
(6  water  quality  PCs  and  31  benthos  PCs  give  186  correlation  coeffi- 
cients) could  lead  to  a number  of  false  correlations  being  accepted 
if  a significance  level  of  0.05  were  adopted.  Furthermore,  it  is 
unclear  to  what  extent  the  benthic  invertebrate  sample  selection 
technique  used  by  Boerger  (1983)  would  affect  the  correlation 
analyses  (see  Section  3.2).  It  was  felt  that,  in  view  of  these  limi- 
tations, it  was  appropriate  to  consider  only  those  correlations 
significant  at  P < 0.001,  recognizing  that  even  these  might  be 
spurious. 

In  May,  the  only  significant  (P  < 0.001)  correlations  were 
between  benthos  PC6  and  water  quality  PC2  and  PC5.  Benthos  PC6  is 
primarily  a measure  of  the  abundance  of  Empididae,  Perlodidae,  and 
Simuliidae  (Table  ^),  and  was  positively  correlated  with  water 
quality  PC2  and  negatively  correlated  with  water  quality  PC5.  Water 
quality  PC2  represents  non-f ilterable  residue,  total  phosphate,  and 
several  metals,  while  PC5  represents  primarily  dissolved 
orthophosphate  (Table  2). 

Water  quality  PCs  2 and  5 also  were  correlated  with  benthos 
PCs  in  June.  Water  quality  PC2  was  negatively  correlated  with  ben- 
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Table  19.  Correlations  between  benthic  invertebrate  principal 

component  scores  and  water  quality  principal  component 
scores  in  1981.  Sample  sizes  are  30  for  May,  60  for  Dune 
and  Duly,  and  30  for  August. 


May; 


Benthos 

Water  Quality 

Correlation 

Principal 

Principal 

Coefficient 

r2 

Component 

Component 

(R) 

Significance 

PCI 

PC6 

0.398 

0.158 

P < 0.05 

PC2 

PC3 

0.383 

0.147 

P < 0.05 

PC5 

PCI 

-0.^65 

0.216 

P < 0.01 

PC2 

-0.^35 

0.189 

P < 0.05 

PC^ 

0.386 

0.149 

P < 0.05 

PC5 

0.506 

0.256 

P < 0.01 

PC6 

PCI 

0.570 

0.325 

P < 0.01 

PC2 

0.573 

0.328 

P < 0.001 

PC4 

-0.513 

0.263 

P < 0.01 

PC5 

-0.657 

0.432 

P < 0.001 

Dune: 

Benthos 

Water  Quality 

Correlation 

Principal 

Principal 

Coefficient 

r2 

Component 

Component 

(R) 

Significance 

PCI 

PC2 

-0.502 

0.252 

P < 0.001 

PC3 

-0.^80 

0.230 

P < 0.001 

PC^ 

0.631 

0.398 

P < 0.001 

PC5 

-0.383 

0.147 

P < 0.01 

PC6 

-0.301 

0.090 

P < 0.05 

PC2 

PC^ 

0.260 

0.068 

P < 0.05 

PC3 

PC5 

-0.^50 

0.202 

P < 0.001 

PC4 

PC^ 

0.if28 

0.183 

P < 0.001 

PC5 

PCI 

0.270 

0.073 

P < 0.05 

PC3 

0.341 

0.116 

P < 0.01 

continued 
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Table  19.  Concluded. 


auly: 


Benthos 

Water  Quality 

Correlation 

Principal 

Principal 

Coefficient 

r2 

Component 

Component 

(R) 

Significance 

PCI 

PCI 

-0.437 

0.191 

P < 0.001 

PC2 

-0.768 

0.590 

P < 0.001 

PC3 

0.455 

0.207 

P < 0.001 

PC4 

0.364 

0.133 

P < 0.01 

PC5 

0.330 

0.109 

P < 0.05 

PC2 

PC3 

-0.585 

0.343 

P < 0.001 

PC3 

PCI 

0.369 

0.136 

P < 0.01 

PC5 

0.316 

0.100 

P < 0.05 

PC5 

PC2 

0.545 

0.297 

P < 0.001 

PC3 

-0.344 

0.118 

P < 0.01 

PC5 

-0.342 

0.117 

P < 0.01 

PC6 

PC6 

0.362 

0.131 

P < 0.01 

PC7 

PC6 

-0.303 

0.092 

P < 0.05 

August: 

Benthos 

Water  Quality 

Correlation 

Principal 

Principal 

Coefficient 

r2 

Component 

Component 

(R) 

Significance 

PC1 

PC4 

0.462 

0.213 

P < 0.05 

PC6 

PC1 

0.431 

0.186 

P < 0.05 

PC2 

0.386 

0.149 

P < 0.05 

PC8 

PC5 

0.394 

0.155 

P < 0.05 
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thos  PC1  and  water  quality  PC5  was  negatively  correlated  with  benthos 
PC3.  Benthos  PC1  in  3une  represents  a relatively  large  number  of 
benthic  taxa,  most  of  which  are  positively  loading,  but  is  primarily 
a measure  of  Chironomidae,  T ricorythodes,  Caenis,  Isonychia,  Hepta- 
qenia,  Cloeon  sp.1,  Oligochaeta,  and  Ophioqomphus  (Table  8).  Benthos 
PC3  represents  Oligochaeta,  Baetis  sp.  A,  and  Simuliidae,  which  are 
positively  loading  taxa,  as  well  as  Ameletus,  Cloeon  sp.  1,  and 
Caenis,  which  have  negative  loadings.  There  were  also  significant 
correlations  between  benthos  PC1  and  water  quality  PCs  3 and  4;  a 
negative  correlation  with  PC3;  and  a positive  correlation  with  PC4. 
Water  quality  PC3  is  a measure  of  potassium,  copper,  zinc,  and  lead, 
while  PC4  represents  primarily  sodium,  chloride,  and  ammonia. 
Benthos  PC4  in  Dune,  which  represents  Ceratopogonidae,  Caenis, 
Brachycercus,  and  Ophiogomphus  (with  positive  loadings)  and  Hydro- 
psychidae,  Ameletus,  and  Heptagenia  (with  negative  loadings),  was 
correlated  positively  with  water  quality  PC4. 

In  3uly,  water  quality  PC1 , which  represents  dissolved 
solids,  was  negatively  correlated  with  benthos  PC1 . Benthos  PC1 
represents  primarily  the  abundance  of  Caenis,  Tricorythodes,  Hepta- 
genia, Ophiogomphus,  Isonychia,  Cladocera,  Cloeon  sp.  1,  and  Collem- 
bola,  all  of  which  are  positively  loading  taxa  (Table  12).  Water 
quality  PC2  (non-f ilterable  residue,  total  phosphate,  and  some 
metals)  was  also  negatively  correlated  with  benthos  PCI,  and  water 
quality  PC3  (potassium,  copper,  zinc,  and  lead)  was  positively 
correlated  with  benthos  PC1 . Water  quality  PC2  was  positively 
correlated  with  benthos  PC5,  and  water  quality  PC3  was  negatively 
correlated  with  benthos  PC2.  Benthos  PC5  in  July  is  a measure  of 
Ostracoda,  Chironomidae,  and  Collembola  (positively  loading)  and 
Baetis  sp.  A and  Isonychia  (negatively  loading).  Benthos  PC2 
represents  primarily  Cloeon  sp.  1 and  Cladocera  (with  positive 
loadings)  and  Isonychia,  Hydropsychidae,  Chironomidae,  Tricorythodes, 
Baetis  sp.  X,  and  Heptagenia  (with  negative  loadings). 

With  some  exceptions,  water  quality  principal  components  1, 
2,  4,  and  5,  which  accounted  for  9 of  the  12  correlations  significant 
at  P < 0.001,  show  consistent  differences  between  left  and  right 
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banks  of  the  Athabasca  River  (Figures  10  to  15).  These  differences 
are  evidence  that  tributary  waters  entering  on  the  right  bank  of  the 
Athabasca,  especially  the  Clearwater  River,  do  not  fully  mix  with 
mainstream  waters  through  much  of  the  study  area.  The  benthos  PCs 
correlated  with  water  quality  PCs  1,  2,  4,  and  5 may  reflect  the 
influence  only  of  drift  and  colonization  of  invertebrates  within  this 
discrete  water  mass,  rather  than  any  direct  influence  of  the  chemical 
constituents  themselves. 

Water  quality  PC3  shows  no  consistent  difference  between 
right  and  left  banks  of  the  Athasbasca  River  (Figure  13).  It  is  not 
clear  what  the  significant  correlations  (P  < 0.001)  between  this 
principal  component  and  three  benthos  PCs  means.  Water  chemistry  PC3 
scores  are  influenced  most  strongly  by  concentrations  of  lead,  potas- 
sium, zinc,  and  copper.  Lead,  zinc,  and  especially  copper,  can  be 
toxic  to  aquatic  life,  so  high  PC3  scores  would  be  expected  to  corre- 
late negatively  with  positively  loading  elements  of  benthic  inverte- 
brate principal  components.  In  at  least  one  case,  however,  there  was 
a positive  correlation  (with  benthos  PCI  in  July).  Therefore,  toxic 
effects  would  seem  to  be  ruled  out.  There  is  no  clear  explanation 
for  the  correlations  with  water  quality  PC3. 
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4.  CONCLUSIONS 

There  are  strong  correlations  apparent  among  several  of  the 
water  quality  parameters  considered  in  the  present  study.  Most  of 
the  major  ions  (calcium,  magnesium,  sodium,  chloride,  sulphate,  and 
bicarbonate)  are  strongly  correlated  with  each  other  and  with 
specific  conductance,  total  alkalinity,  and  filterable  residue. 
Potassium  concentrations,  however,  appear  to  be  independent  of  the 
concentrations  of  other  major  ions.  There  are  also  components  of 
sodium  and  chloride,  attributable  to  the  Clearwater  River  and  other 
east  bank  tributaries,  that  are  not  strongly  related  to  the  concen- 
trations of  other  ions.  All  of  the  metals  except  lead  and  mercury 
are  strongly  associated  with  non-f ilterable  residue  and  total  phos- 
phate. There  is  also  some  suggestion,  however,  that  there  may  be 
components  of  zinc  and  copper  that  are  not  associated  with  non- 
filterable  residue. 

With  respect  to  those  water  quality  parameters  included  in 
the  present  study,  there  is  no  evidence  that  effluent  from  the  Suncor 
plant  has  a large  or  consistent  effect  on  the  water  quality  of  the 
Athabasca  River  downstream  from  the  development.  The  suggestion  is 
that  any  effects  on  water  quality  are  short-term  in  nature  and  have 
not  resulted  in  changes  that  persist  for  long  periods  of  time. 

The  major  differences  in  water  quality  within  the  study 
area  are  between  the  left  and  right  sides  of  the  Athabasca  River 
downstream  from  Fort  McMurray.  These  differences  are  attributable 
primarily  to  the  influences  of  the  Clearwater  River  and  other  right 
bank  tributaries,  and  involve  mainly  sodium  and  chloride  concentra- 
tions. 

There  are  some  differences  in  benthic  invertebrate 
abundance  and  community  composition  within  the  study  area  that  appear 
to  be  related  to  nutrient  enrichment  due  to  Fort  McMurray  sewage  eff- 
luent (see  Section  3.2.5.1).  There  is  no  evidence  of  large  differ- 
ences in  benthic  invertebrate  populations  between  stations 
immediately  upstream  and  downstream  of  the  Suncor  development  that 
can  be  attributed  to  the  Suncor  effluent.  Most  of  the  observed 
differences  in  benthic  invertebrate  abundance  or  community  composi- 
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tion  are  differences  between  the  left  and  right  sides  of  the 
Athabasca  River,  and  are  probably  due  to  the  influences  of  right  bank 
tributaries  and,  at  least  for  several  kilometres  downstream  from  Fort 
McMurray,  to  the  effects  of  effluent  from  the  sewage  treatment  plant. 

The  abundance  of  a variety  of  benthic  invertebrate  taxa 
appears  to  be  correlated  with  several  water  quality  parameters.  In 
most  cases,  these  correlations  are  with  water  quality  parameters  that 
show  differences  between  the  left  and  right  sides  of  the  Athabasca 
River.  The  observed  correlations  therefore  may  be  due  to  drifting 
of,  and  colonization  by,  invertebrates  from  the  Clearwater  River,  and 
possibly,  other  right  bank  tributaries,  rather  than  to  any  direct 
influence  of  water  quality. 

The  results  of  analysis  of  benthic  invertebrate  data 
presented  in  this  report  should  be  considered  tentative  or  incon- 
clusive. The  method  used  for  selecting  benthic  samples  adversely 
affects  the  validity  of  statistical  analyses  (see  Section  3.2),  and 
the  reliability  of  the  results  is  therefore  questionable. 
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5.  RECOMMENDATIONS 

Within  the  group  of  29  water  quality  parameters  examined  in 
detail  in  the  present  study,  there  would  appear  to  be  few  opportuni- 
ties to  limit  the  number  of  parameters  measured  to  a smaller  group 
and  still  adequately  characterize  water  quality.  The  ions  of 
calcium,  magnesium,  sulphate,  and  bicarbonate  are  all  strongly  corre- 
lated with  each  other  and  with  specific  conductance,  total  alkalin- 
ity, and  filterable  residue.  It  may  not  be  necessary,  therefore,  to 
obtain  measurements  on  all  of  these  parameters.  In  particular, 
specific  conductance,  total  alkalinity,  and  filterable  residue  would 
appear  to  provide  little  information  in  addition  to  that  already 
provided  by  measurements  of  other  parameters. 

The  pattern  of  correlations  observed  for  measurements  of 
nitrogen,  phosphate,  and  organic  carbon  suggest  that  continued 
measurement  of  all  these  would  be  appropriate.  Similarly,  phenolic 
material  and  oil  and  grease  are  largely  independent  of  other  parame- 
ters and  of  each  other,  and  they  also  should  be  measured  in  future 
water  quality  monitoring. 

Although  there  are  strong  correlations  among  most  of  the 
metals  (except  lead  and  mercury),  the  results  of  two  principal 
component  analyses  (see  sections  3.1.1  and  3.1.2)  are  not  entirely 
consistent  in  this  regard.  It  would  not  appear  practical,  at  the 
present  time,  to  use  measurements  of  some  of  the  metals  as  indicators 
of  the  probable  levels  of  other  metals. 

The  number  of  water  quality  monitoring  stations  and  their 
locations  in  the  study  area  would  appear,  for  the  most  part,  to  be 
adequate  for  the  present  level  of  development  in  the  area.  One 
additional  station,  downstream  from  Fort  McMurray  but  closer  to  the 
city,  probably  would  be  useful.  Sampling  at  monthly  intervals,  which 
has  usually  been  the  schedule  used,  appears  to  be  generally  adequate 
if  only  continuous  effluent  discharges  are  involved,  but  more 
frequent  sampling  during  periods  of  rapid  change  would  seem 
desirable.  If  the  application  of  multivariate  statistical  methods  is 
anticipated,  it  is  important  that  all  water  quality  parameters  are 
determined  for  all  samples. 
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When  benthic  invertebrates  are  sampled  for  the  purposes  of 
a monitoring  study,  and  when  the  effluent  discharges  of  interest  are 
continuous,  three  sampling  periods  during  the  year  (early  spring, 
mid-summer,  and  fall)  should  be  sufficient.  Many  closely  spaced 
sampling  intervals  are  not  particularly  useful  unless  a precise 
description  of  seasonal  changes  is  required.  In  general,  it  would  be 
better  to  reduce  the  number  of  sampling  periods  and  increase  the 
sample  sizes.  Sample  sizes  of  three,  although  common  in  benthic 
invertebrate  studies,  are  minimal,  and  it  is  unlikely  that  less  than 
catastrophic  changes  in  benthic  populations  would  be  detectable  with 
such  a small  number  of  replicate  samples.  It  is  essential,  for  the 
valid  application  of  statistical  methods,  that  random  sampling 
procedures  be  employed  when  collecting  benthic  invertebrate  samples. 

For  the  purposes  of  examining  relationships  between  benthic 
invertebrates  and  water  quality,  it  would  be  desirable  to  conduct  a 
study  at  a location  where  a range  of  water  quality  conditions  exists 
and  where  possible  relationships  are  not  confounded  to  a great  extent 
by  the  effects  of  drift  and  colonization  of  benthos  from  major  tribu- 
taries to  the  stream  being  studied.  Benthic  populations  in  tributary 
streams  may  be  quite  different  for  a variety  of  reasons  unrelated  to 
water  quality.  In  addition,  it  is  clear  that  where  benthic  inverte- 
brate populations  may  be  influenced  by  water  quality,  the  influence 
takes  place  over  a period  of  time  and  is  not  constant  during  this 
time.  Water  quality  measurements  taken  at  the  same  time  as  benthic 
samples  may  not  be  particularly  relevant  in  terms  of  the  observed 
benthic  populations.  An  alternative  approach,  that  might  be  more 
informative,  would  be  to  integrate  water  quality  measurements  over 
some  period  of  time  prior  to  sampling  benthic  invertebrates.  These 
integrated  values  would  then  be  used  as  an  indication  of  the  water 
quality  conditions  to  which  a particular  site  had  been  exposed 
previously.  This  approach  would  require  regular  monitoring  of  water 
quality  for  a period  of  time  prior  to  the  sampling  of  benthic 
invertebrates.  It  is  suggested  that,  at  least  as  a first  guess,  a 
time  period  of  one  full  year  in  advance  of  benthic  sampling  may  be 
appropriate. 
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7.  APPENDICES 

7.1  WATER  QUALITY  DATA,  1976  to  1983 

Figures  16  to  83  include  plots  of  68  water  quality 
parameters  measured  at  seven  stations  on  the  Athabasca  River  for  the 
years  1976  through  1983.  The  determination  of  location  in  the  river 
channel  was  based  on  distances  to  the  left  and  right  banks  (looking 
downstream)  as  recorded  in  the  NAQUADAT  data  files.  Samples  located 
anywhere  within  the  left  34%  of  the  river  width  were  considered  left 
bank  samples,  those  within  the  right  34%  were  considered  right  bank 
samples,  and  those  within  the  central  32%  were  considered  mid-channel 
samples.  For  105  of  the  538  water  samples,  there  was  insufficient 
information  in  the  NAQUADAT  files  to  determine  sampling  location. 
These  samples  were  assumed  to  be  mid-channel  samples.  Curves  on  the 
graphs  were  interrupted  whenever  the  time  period  between  samples  was 
greater  than  four  months. 
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Figure  16.  Water  temperature  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1983. 
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Figure  17.  Field-measured  pH  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1983. 
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Figure  18.  Laboratory-measured  pH  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1963. 
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Figure  19.  Dissolved  oxygen  at  seven  stations  on  the  Athabasca  River 
from  1976  to  1983. 


continued 


mg/L  mg/L  mg/L 


105 


15.0  n 


10.0- 


5.0- 


STnTION  DR0206 


0.0 

15.0  n 


10.0-  / 


5.0- 


0.0 


I 

I 

I 


H ■ I ■ 1 1 ■ I ■ ■ ■ ■ V r.i . ■ f 1 r I-.  |I  ■ ITTI  ■ ri  np 

STRTION  DR02Q7 


\( 


r- 


\J 


> 

V' 


Parameter:  08102P 
OlseoLved  Oxggen 

Fiaure  19.  Concluded. 


ujo/gri  ujo/grl  uio/grl 


106 


1976  1977  1978  1979  1980  1981  1982  1983 


LEGEND 

Left  bank 

Parameter:  0204 IF  rLqhi  bank 

Specific  Conductance  mid-channel 


Figure  20.  Field-measured  specific  conductance  at  seven  stations  on 
the  Athabasca  River  from  1976  to  1983. 
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Figure  21 . Laboratory-measured  specific  conductance  at  seven 
stations  on  the  Athabasca  River  from  1976  to  1983. 
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Figure  22.  Calcium  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  23.  Kiagnesium  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  24.  Total  hardness  at  seven  stations  on  the  Athabasca  River 
from  1976  to  1983. 
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Figure  25.  Potassium  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  26.  Sodium  at  seven  stations  on  the  Athabasca  River  from  1976 
to  1983. 
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Figure  27.  Chloride  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  28.  Sulphate  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  29.  Silica  at  seven  stations  on  the  Athabasca  River  from  1976 
to  1983. 
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Figure  30.  Fluoride  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  31.  Total  alkalinity  at  seven  stations  on  the  Athabasca  River 
from  1976  to  1983. 
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Figure  32.  Phenolphthalein  alkalinity  at  seven  stations  on  the 
Athabasca  River  from  1976  to  1983. 
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Figure  32.  Concluded. 
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Figure  33.  Bicarbonate  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  34.  Carbonate  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  35.  Total  inorganic  carbon  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1983. 
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Figure  36.  Free  carbon  dioxide  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1983. 
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Figure  37.  Filterable  residue  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1983. 
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Figure  38.  Filterable  residue  (fixed)  at  seven  stations  on  the 
Athabasca  River  from  1976  to  1983. 
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Figure  39.  Non- filterable  residue  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1983. 
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Figure  40.  Non-filterable  residue  (fixed)  at  seven  stations  on  the 
Athabasca  River  from  1976  to  1963. 
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Figure  40.  Concluded. 
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Figure  41.  Turbidity  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  ^2.  Chlorophyll  a at  seven  stations  on  the  Athabasca  River 
from  1976  to  1983. 
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Figure  43.  Total  organic  carbon  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1983. 
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Figure  44.  Dissolved  organic  carbon  at  seven  stations  on  the 
Athabasca  River  from  1976  to  1983. 
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Figure  45.  Tannin  and  lignin  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1983. 
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Figure  46.  Humic  acid  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  46.  Concluded 
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Figure  4-7.  Fulvic  acid  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  48.  Water  colour  at  seven  stations  on  the  Athabasca  River 
from  1976  to  1983. 
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Figure  49.  Nitrate  plus  nitrite  nitrogen  at  seven  stations  on  the 
Athabasca  River  from  1976  to  1983. 
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Figure  50.  Nitrite  nitrogen  at  seven  stations  on  the  Athabasca  River 
from  1976  to  1983. 
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Figure  50.  Concluded. 
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Figure  51 . Nitrate  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  52.  Ammonia  nitrogen  at  seven  stations  on  the  Athabasca  River 
from  1976  to  1983. 
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Figure  53.  Total  Kjeldahl  nitrogen  at  seven  stations  on  the 
Athabasca  River  from  1976  to  1983. 
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Figure  54.  Total  phosphate  phosphorus  at  seven  stations  on  the 
Athabasca  River  from  1976  to  1983. 
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Figure  55.  Orthophosphate  phosphorus  at  seven  stations  on  the 
Athabasca  River  from  1976  to  1983. 
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Figure  56.  Surfactants  at  seven  stations  on  the  Athabasca  River  from 
from  1976  to  1983. 
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Figure  57.  Phenolic  material  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1983. 
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Figure  58.  Oil  and  grease  (determined  by  petroleum  ether  extraction) 
at  seven  stations  on  the  Athabasca  River  from  1976  to 
1983. 
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Figure  59.  Oil  and  grease  (determined  by  freon  extraction)  at  seven 
stations  on  the  Athabasca  River  from  1976  to  1983. 
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Figure  60.  Hydrocarbons  (alkanes)  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1983. 
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Figure  61.  Dissolved  sulphide  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1983. 
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Figure  62.  Cyanide  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  63.  Iron  at  seven  stations  on  the  Athabasca  River  from  1976 
to  1983. 
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Figure  64.  Aluminum  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  65.  Manganese  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  66.  Boron  at  seven  stations  on  the  Athabasca  River  from  1976 
to  1981. 
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Figure  67.  Copper  at  seven  stations  on  the  Athabasca  River  from  1976 
to  1983. 
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Figure  68.  Zinc  at  seven  stations  on  the  Athabasca  River  from  1976 
to  1983. 
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Figure  69.  Lead  at  seven  stations  on  the  Athabasca  River  from  1976 
to  1983. 
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Figure  69.  Concluded 
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Figure  70.  Mercury  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  71.  Arsenic  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  71.  Concluded. 
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Figure  72.  Selenium  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  73.  Vanadium  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  74.  Chromium  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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Figure  76.  Cobalt  at  seven  stations  on  the  Athabasca  River  from  1976 
to  1983. 
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Figure  76.  Concluded. 
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Figure  77.  Cadmium  at  seven  stations  on  the  Athabasca  River  from 

1976  to  1983. 
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Figure  77.  Concluded. 
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Figure  78.  Silver  at  seven  stations  on  the  Athabasca  River  from  1976 
to  1983. 
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Figure  79.  Fecal  coliform  bacteria  at  seven  stations  on  the 
Athabasca  River  from  1976  to  1983. 
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l-igure  79.  Concluded. 
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Figure  80.  Total  coliform  bacteria  at  seven  stations  on  the 
Athabasca  River  from  1976  to  1983. 
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Figure  81 . Standard  bacterial  plate  counts  at  seven  stations  on  the 
Athabasca  River  from  1976  to  1983. 
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Figure  82.  Chemical  oxygen  demand  at  seven  stations  on  the  Athabasca 
River  from  1976  to  1983. 
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Figure  82.  Concluded 
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Figure  83.  Beryllium  at  seven  stations  on  the  Athabasca  River  from 
1976  to  1983. 
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7.2  BENTHIC  INVERTEBRATE  ABUNDANCE  DATA 

Benthic  invertebrate  abundance  data  for  16  sites  on  the 
Athabasca  River  during  the  spring  and  summer  of  1981  are  presented  in 
tables  20  to  23.  The  data  are  from  Boerger  (1983),  and  for  each  of 
May,  Dune,  and  Ouly,  the  samples  from  two  sampling  periods  were 
pooled  for  calculating  means  and  standard  deviations.  The  August 
data  includes  samples  from  only  one  sampling  period. 


Table  20,  Benthic  Invertebrate  abundance  (number/m^)  at  16  sites  on  the  Athabasca  River  between  May  13 
and  May  29,  1981.  Values  are  means  and  standard  deviations  based  on  6 samples  from  each 
station. 
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Table  21.  Benthic  Invertebrate  abundance  ( number/ sites  on  the  Athabasca  River  between  Oune  9 
and  June  2^,  1981.  Values  are  means  and  standard  deviations  based  on  6 samples  from  each 
station. 
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and  July  22,  1981.  Values  are  means  and  standard  deviations  based  on  6 samples  from  each 
station. 
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Table  23.  Benthic  Invertebrate  abundance  (number/m2)  at  16  sites  on  the  Athabasca  River  between  August 
18  and  August  19,  1981.  Values  are  means  and  standard  deviations  based  on  3 samples  from 
each  station. 
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